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Objective 


To develop an improved model that can be incorporated into the FSCBG, AGDISP, or any other spray model to 
predict changes in droplet size and concentration due to evaporation. 


Background 


A critical factor in the aerial spray predictions models is the evaporation rate associated with the spray droplets. 
Evaporation affects the droplet size and hence the terminal velocity of the drop. Consequently, evaporation has 
significant impact on the time and space distribution of the spray. A study by Dr. Banaugh, in 1979, revealed that 
most of the theoretical research was done for pure liquid drops. In order to have information for the spray models, a 
study was sponsored at Colorado State University to measure evaporation rate of droplets from 13 pesticide 
mixtures under simulate free fall conditions. No further use was made of the windtunnel and it was removed and 
dismantled. Both the AGDISP and FSCBG are configured to accept data from the simulated tests or to calculate 
evaporation rate for pure liquids when certain physical properties are known. 


This year a new study was sponsored to review and update the status of droplet evaporation research and 
development. Eight organizations have been identified that have done recent testing of droplet evaporation. The 
U.S. Army, CRDEC, at Aberdeen Proving Ground has developed the most complete and thorough evaporation 
model for multi-component mixtures, including emulsions and suspensions. NUSSE4 (the CRDEC model) can 
simulate up to five components. The five components can include two volatile and three non-volatile components. 
Only limited validation of the model predictions have been made to date. 


Since the pesticide tank mixes have changed, we no longer have valid experimental results for the spray model, nor 
do we have multi-component models operational. 
Methods 


We plan to obtain the NUSSE3 and NUSSE4 models and make them operational on a Forest Service computer. 
We will determine the physical constants needed to exercise the model and obtain those constants from the 
pesticide manufacturers if possible. We plan to obtain experimental measurements of new tank mixes and compare 
those results to the model predictions. Comparisons of previous experimental results will be made to existing 
models. Revisions to the model will be proposed based on comparisons to experimental measurements. 


Revision 


Early examination of the NUSSE models revealed that their current use is limited to certain specialized compounds 
of interest to the Army. Further investigation of the single component evaporation model, already a part of the 
FSCBG code required measurements of physical properties that are difficult to measure even in a laboratory. 


We are therefore recommending that we pursue an alternative approach. This is, contract for measurement of 
evaporation of single droplets. This is the method that was developed under contract by Colorado State University 
in the early 1980's. (Appendix A). The direct measurement of a matrix of droplet diameters, temperatures, relative 
humidities, and tank mixes can be input directly into either AGDISP or FSCBG. 


The system developed by Colorado State University is documented, but has been dismantled. We have not been 
successful in identifying another laboratory that has an operational system to measure droplet evaporation directly. 


Most of the funding for this project is intact and can be applied to a contract to develop a system and conduct 
measurements of evaporation rates of single droplets. 





Workplan Revised July 1, 1990 
Selection and Verification of Spray Droplet Evaporation Model 











Assignment 1990 1991 
1/2/3) 4 Tec 
Compare AGDISP and FSCBG Predictions Thompson X| X 
to 1984 Laboratory Measurements Teske 
Measure Evaporation Rates of Single Droplets Contractor XHEX X| X 
Ekblad 
Barry 
Write request for proposals 
Evaluate proposals 
Award contract 




















Monitor contract 











Prepare Quadratic Coefficients Contractor x 
for AGDISP and FSCBG inputs 














Appendix—Determination of Evaportion Rates of Pesticide Droplets 










oe 
ect a 
oar Pa. 
: une 1 


cm 
pare AGOISP 
~ re 
yap 
Jy ioe ree 
aaa 
vi ra 
: . + 
“at : 
or 
7H a 


> a] ™ 
al 
aft. 
+R 
— ape 
ara: b 
Ate, ve Wig 
f- 
: 
’ @ 
: 


rte watts. 

pear . 

Ment i 7 
ribet | 


Wvangre Quis tests 
; 
(AGI | 
> & 
Sao 
oe 
SP ad 
i tah 


Determination of 
Evaporation Rates of 
Pesticide Droplets 


FINAL REPORT 


Prepared by 

Rick S. Dennison 

Civil Engineering Department 
Colorado State University 


James B. Wedding, Ph.D. 
Director, Aerosol Science Laboratory 
Affiliated with the 

Research Institute of Colorado 


Prepared for 

U. S. Department of Agriculture 
Forest Service 

Equipment Development Center 
Missoula, MT 59801 


Robert B. Ekblad 


Project Leader 


Contract No. 53-0343-8-105 
(Work under this contract was completed in July 1982) 


Pesticide Precautionary Statement 


This publication reports research involving pesticides. It does 
not contain recommendations for their use, nor does it imply 
that the uses discussed here have been registered. All uses of 
pesticides must be registered by appropriate State and/or 
Federal agencies before they can be recommended. 


CAUTION: Pesticides can be injurious to humans, domestic 
animals, desirable plants, and fish or other wildlife-if they 
are not handled or applied properly. Use all pesticides selec- 
tively and carefully. Follow recommended practices for the 
disposal of surplus pesticides and pesticide containers. 
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FOREWORD 


This report is published as a part of a USDA Forest Service program to 
improve aerial application of pesticides, specifically by using pesticides 
and delivery systems tailored to the forest environment. The program is 
conducted jointly by the Equipment Development Center, Missoula, Mont., 
and Forest Pest Management Staff, Washington Office at Davis, Calif., 
under the sponsorship of State and Private Forestry. 


Details of the aerial application improvement program are explained in two 


Forest Service reports, A Problem Analysis: Forest and Range Aerial Pesticide 
Application Technology (Equipment Development Center Rpt. 7934 2804, July 1979, 


Missoula, Mont.) and Recommended Development Plan for An Aerial Spray Planning 
and Analysis System (Forest Pest Management Rpt. FPM 82-2, February 1982, 
Davis, §Caliis)% 


The system of computer models developed to optimize spray program design 
and operation and evaluate assessment of environmental risk posed by aerial 
Spray operations must include predictions of evaporation rates of spray 
droplets. The classic theoretical evaporation equations are suitable for 
pure liquids, such as water, but have not been extended to more complex 
mixtures, emulsions, and suspensions ordinarily used in spraying. The 
alternate to the use of theoretical equations is experimental measurement 
under controlled conditions. 


This report contains the results of an extensive series of experimental 
measurements of droplet evaporation rates. Fourteen tank mixes representing 
a variety of pesticides, a commonly used antidrift additive, and water were 
examined. 


Four sizes of droplets were examined under simulated free fall conditions 
for a wide variety of temperature and relative humidity conditions. Enough 
data were taken to insure that a reasonably accurate estimate of other 
droplets and conditions could be made by interpolation. Some examples of 
evaporation rates are included in the report but the majority of the data 
are stored on magnetic tape at the USDA Fort Collins Computer Center. 
Complete instructions for retrieving the data are included in the report. 


The study and measurements were conducted by the Aerosol Science Laboratory 
of the Engineering Research Center at Colorado State University under 
contract to the USDA Forest Service. 


The droplet evaporation facility has been retained at the Aerosol Science 
Laboratory and is available for analysis of other tank mixes on a contract 


basis. 


A report on Equtpment Development and Test Project 2664, Optimum-Stze Spray 
Droplets, funded by the Forest Pest Management Staff. 
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Final Report 
DETERMINATION OF EVAPORATION RATES OF 
PESTICIDE DROPLETS 

An experimental study was conducted to measure the rates of mass 
transfer from 13 pesticide mixtures plus water under simulated free 
fall conditions. Droplets were captured on a wire and air that was 
conditioned at various temperatures and relative humidities was passed 
by the pendant droplet in a closed loop wind tunnel facility. Droplets 
having nominal diameters of 50, 100, 250, and 400 microns generated 
uSing a vibrating reed aerosol generator were sized at each time step 
and the wind tunnel velocity adjusted to the corresponding terminal 
velocity. The droplet diameter was determined using a Reticon solid 
State camera and an on-line Texas Instruments microcomputer. 

The data for Water were: used to obtain, an equation of the 
Sherwood number (Sh) as a function of the Reynolds (Re) and Schmidt (Sc) 
numbers. The resulting correlation equation, 
ois) lesson eee Cy ee ain agreed well with previous authors. The 
data obtained for the other test Solutions were reduced to plots of percent 
of mass loss versus time. These curves allowed direct application of 
the results obtained in this study. 

The measurement and control of the air humidity presented a 
possible source of error. However, the data reported generally repre- 
sent the mass transfer of freely falling spray droplets. An overall 
error of approximately 9% resulted from this experimental procedure 
anda tact Lieys 

Rick S. Dennison 

Civil Engineering Department 
Colorado State University 
Fort. Collins, Colorado 80523 


spring, 1982 
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Chapter 1 


INTRODUCTION 


The control and management of various forest insects has been accom- 
plished recently with the aid of aerially disseminated pesticides. In 
order to ontimize this use of the vesticide spray, several parameters must 
be determined. These parameters include initial droplet size distribution, 
volatility of the material, height and width of the spray release, and the 
cloud behavior affected by the aircraft and atmospheric fluid dynamics. 
Complete knowledge of the above could allow the efficient application of 
pesticides from aircraft. The purpose of the present study was to examine 
several spray solutions and determine their corresponding time rate evap- 
oration histories. Droplets were to be examined under simulated free fall 
conditions for a variety of temperature and relative humidity combinations 


of the passing airstream. 


Chapter 2 
THEORETICAL CONSIDERATIONS 
2.1 Evaporation from Stationary Droplets 
Early work with mass transfer was performed by Maxwell (8) for 
the evaporation of a spherical droplet. In this study, Maxwell devel- 


oped an expression; for. the,ratesofemass lossfofear droplet 


dm Cece 
(Ge), =2247nDr ar G24) 


Equation 2.1 assumes the radii of the droplets are large compared, toe,the 
mean free path of the surrounding fluid. This expression definés the 
droplet evaporation for the simplest case, that of a quiescent, homo- 
geneous medium. Langmuir (7) also derived equation 2.1 in an independent 
study performed years later. “Integration of, equation 2 41° and sche 
substitution ‘of boundary conditions yields: 

(ie), = -4raD(Cp co (2.2) 


whereas the droplet radius, Ce 1s Chesconcentration of vapor at, the 


0 


surface, and C| is the free stream concentration. Assuming the vapor 
of the evaporating liquid obeys the ideal gas laws, the concentration of 
the vapor can be expressed as a partial vapor pressure, P, such that: 


Cat (28S) 
RT 


Substitution of equation 2.5 into equation 2.2 reduces the integrated 
Maxwell equation: 


(=) _ _ 4naDM, & . =| 
dt 0 or To i (2.4) 





where P) and T) describe the conditions at the droplet surface and P. 


and T_ correspond to the surrounding medium. 
2.2 Effects of Forced Convection ; 

The mass transfer process of a droplet in a motionless, uniform 
medium is completely defined by equation 2.4. When a relative motion 
between the droplet and the medium is introduced, the effects of forced 
convection are described by the Sherwood number, Sh. This dimensionless 


chemical engineering parameter is defined as: 


(i) 
Sie Wits eke (2.5) 
DS (C,-C.) 
(aE) 
where \dt : is the convective mass rate and S is the droplet surface 


area. Equation 2.5 can be rearranged to include the integrated form of 
Maxwell's equation (2.2) for motionless evaporation: 
($2) (St) 
2a\dt dt as 
< — fae) 


. [-4na“D(C,)-C,)] : 2 


Sh 
(ae 


0 
Similarly, physical properties of heat and mass transfer allow the effects 


of forced convection on the rate of heat transfer to be described by the 


Nusselt number: 


Goes) 


Close examination of equations 2.6 and 2.7 reveals that Sh and Nu are 
defined as wind factors. One half of these values are multiplication 
factors used with the motionless rates of transfer to obtain the larger 
rates experiences with forced convection. When a relative motion 
between the droplet and the medium does not exist, Sh and Nu are reduced 
to a constant value of two. This value defines the theoretical lower 


Limiteort ohv and Nu: 


Chapter 3 
Review of Literature 

The effects of forced convection on the heat and mass transfer of 
spherical droplets has been examined at great extent. Theoretical and 
experimental solutions for the Sherwood and Nusselt numbers have been 
produced for a wide range of ventilation velocities and droplet diameters. 
The following review will consider solutions dealing with mass transfer 
only. ““Equations{similarvto these “describing ‘the waneeee of heat are 
often presented by authors due to the analogous physical procedures 
involved. 

Limited success has been obtained using strict theoretical solutions 
to describe the rates of transfer. This limitation is caused by the 
difficulty in solving the Navier-Stokes equation of fluid flow around 
the droplet simultaneously with the differential equations of heat and 
mass transfer. However, several authors have obtained solutions by 
making assumptions about the physical process. Froessling (4), as cited 
in Fuchs (S), used a boundary layer theory to obtain equations for a 
blunt-nosed body of revolution. With this solution, he determined that 
the Sherwood number could be expressed as a function of two additional 
dimensionless parameters, the Reynolds (Re) and Schmidt (Sc) porbaoh 


These dimensionless groupings are defined as follows: 


Iq VP 
2aV"a ipa (3eN) 


m1 p ,D 


: : 172 
Specifically, Froessling presented Sh as being proportional to Re ae 


Subsequent experiments using droplets of water, aniline, nitrobenzene, 
and naphthalene were carried out by Froessling. These experiments were 
performed with constant ventilation velocities using droplets suspended 
from glass fibers. An environmental temperature of 20°C was used in 
everyutest. “dntall ‘cases; thecdatagwerencorrelat edaby.: 


She 2. Oat BRei/? ge) ® 2) 


where 8, called the mass transfer coefficient, was found to be 0.552. 
This equation remained consistent for a wide range of Re, from 2.3 to 
1280. 

Expressions duplicating equation 3.2 have been presented by inde- 
pendent studies. Fuchs, using a similarity principle, determined that 
the Sherwood number averaged over the entire droplet surface was pro- 


portional to Rea ee 


He combined this result with the limiting 
case of She=.2 wheneRepg=- 0) to, obtaingequat iony3.Z.geThe author, defines 
B aS a constant to be determined experimentally, but does not include 
data to calculate a value. Ranz and Marshall (10) also obtained an 
expressiongidenticaljto equation!3:2 JeThesauthors used Froessling's 
boundary layer equations in combination with the equations of heat and 
mass transfer to obtain a set of dimensionless partial differential 
equations. From these equations, a theoretical and functional form for 
the Sherwood number was obtained, Sh = Sh(Re, Sc). Experimental data 
for droplets of water and benzene obtained by Ranz and Marshall were 
correlated with equation Oe The experimental procedure consisted of 
measuring the rate of feed necessary to maintain a constant diameter 
droplet suspended on a microburet. Dry air with temperatures up to 


O ‘ 
220 C was used as the medium. The mass transfer coefficient, 8B 


determined by this data was equal to 0.6. 


Kinzer and Gunn (6) developed a theory that considered the transient 
transfer of vapor to successive packets of fresh environmental air as 
the droplet moved relative to the medium. In this theory, the authors 
sum the transfer contributions of all the packets that contact the 
droplet. In doing so, estimations of the contact area, A, and the time 
of contact for each packet, t, were required. The vapor transferred to 


each fresh packet was given by: 


At (S) ; (3.3) 


Wiere™x” PS “Wire Vapor density Patio. and ($5) is’ the "sradient at ‘tHe 
a 


droplet surface. The values obtained for each packet by equation 3.3 had 
to be summed over a time interval, ty = 2a/v, in order to determine the 
total tran: fer. The*final form of ‘this theoretical approach was given 


ase 


Shieeaead (ls 0 t+ oP (Re Sc/4n) *] (3.4) 
The value of F, called the "ventilation factor'', was defined as the 
ratio of the sum of the terms of equation 3.3 for the ventilated drop- 
ie Cero i nessame duantityelorea droplet vat rest. (The ventilation factor 
depended on the ventilation velocity and had to be determined by direct 
observation. Experiments were performed by Kinzer and Gunn with water 
droplets using a wide range of temperatures (0 - 40°C) and humidities 
(10 - 100%) under i fall conditions. Three methods of size measure- 
ment were used in these experiments for different ranges of droplet 
diameters. The method of successive photographs was used to record 
progressive positions of”droplets ranging from 10 to 140 microns in 
diameter. The droplets were generated at the top end of a well-insulated 


shatt and allowed to fali freely to the bottom.. The terminal settling 


velocities were determined using the photographs and the time between 
each photograph. The droplet diameters were determined from the depen- 
dence of diameter to terminal velocity. Similarly, droplets with 
diameters from 40 to 1000 microns were examined as they fell through a 
series of inducing rings. Droplets were charged as they were produced 
and were traced as they fell through the inducing rings. The charged 
droplets produced an impulse as they passed through each ring. The time 
between each impulse was used to find the terminal velocities and in 
turn the droplet diameters. The third method measured droplet diameters 
from 0.1 to 0.3 cm as a function of height maintained inside a conical 
tube experiencing a constant velocity. The evaporation rates of the 
droplets were obtained by measuring the decreasing heights of the droplets 
at successive time intervals. Comparison of the heights with calibrated 
data produced diameter versus time data. The results of these experi- 
ments were evaluated by equation 3.4 so as to determine values of F for 
Cachedropket diameter, s ne taumiorsspresen leds tiesresulLGseassaaplOtsor 


E/2 


Feversus, Re It should be noted that the results given by Kinzer and 


Gunn indicate that sh-isenotiuaslinear tunction or Ree 


Thiss 1s2a 
complete contradiction to the results given by Froessling and Ranz and 
Marshall correlated by equation 3.2. 

Davies (3) reviewed the studies dealing with the transfer of mass 
from droplets. He presents a step by step procedure used to calculate 
rates of evaporation in still air and for forced convection. In this 
procedure, equation 3.2 was used to describe the effects of forced 


convection. Davies also presented experimental data performed at an 


air temperature of 20% using droplets of furfural and water. These 


droplets were suspended on a fine glass spill in a constant air stream. 
The author used equation 3.2 with 8 = 0.552 to compare with the results 
of the experiments. Errors up to 13% were obtained in this comparison. 
Beard and Pruppacher (1) developed a unique experimental procedure 
to determine the mass transfer from droplets under free fall conditions. 
Experiments with water droplets were performed using a limited range of 
temperatures (22 - 24°C) and humidities (27 - 65%). Droplets ranging 
in diameter from 100 - 375 microns were held stationary by a vertical 
tunnel air stream. The velocities of the air stream required to main- 
tain the stationary position of the droplet were recorded over time. 
The velocity record was then used to determine the rate of change of 
droplet size. The data were correlated with an equation similar to 
equation 3.2. Two distinct linear curves were found over the range of 


pel/ 2 col/> tested. For 1.43 = Re pl < 12.0, the authors found 


that: 
Sh = 1.56 + 0.616 Ree ae (355) 
12 L/S : 
Values of Re SV below 1.43 were described by the relation: 
Sh = 2.0 + 0.216 ret /? eo (326) 


Equation 3.6 was obtained to satisfy the theoretical lower limit, 


Sh = 2.0, along with the experimental data. 


Chapter 4 
SCOPE OF STUDY 

The present study was to determine time rate evaporation histories 
for fourteen solutions. These solutions include thirteen mixtures of 
pesticides or pesticide components used today and water as a reference. 
The Evaporation Rate Facility of the Aerosol Science Laboratory at 
Colorado State University utilized a microscope coupled to a solid state 
camera to measure the droplet diameter directly at each present time 
step. A microcomputer programmed by Aerosol Science Laboratory 
personnel calculated the droplet diameter from the image of the pendant 
droplet as seen by the camera. The microcomputer then calculated the 
terminal settling velocity and adjusted the wind tunnel velocity to 
equal this terminal velocity. The droplet was sized and velocity 
adjusted continuously at each. second of the test. 

The test condition variables for each of the 14 solutions 
included environmental temperature, relative humidity of the air, and 
nominal initial droplet diameter. The testing values used in this study 
for these variables can be seen in Table 1. By changing these conditions, 
the relative effects on the mass transfer from the droplets was 
ascertained. 

Varying the initial droplet diameters allowed the mass transfer 
process of each solution to be described for different sized droplets 


falling at terminal velocities. The wide range of test parameters used 
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Tabierel 


Experimental Testing Conditions 


Environmental Temperatures 
eG 
10°C 
LS 


20°C 


Environmental Relative Humidities 


Initial Droplet Diameters 
50 microns 
100 microns 
250 microns 


400 microns 


1 


in the study allowed most common environmental conditions to be 
matched. Thus, the results produced by this study could be applied 
directly to a proposed spray operation. For a given pesticide, nozzle 
spray droplet distribution, environmental temperature and relative 
humidity, and spray height the mass of pesticide that reaches the 
intended target could be approximated, if one neglects the added 
effects of mixing from the aircraft wake and atmospheric turbulence. 
This would assist in the determination of the effect of the pesticide 
on the target insect population. The total mass that reaches the 
intended target will increase or decrease with variations of the envir- 
onmental conditions. Optimum conditions can then be determined that 
allow the maximum effect of the pesticide on the insect population. A 
cost and effect evaluation could also be applied to optimize the entire 
procedure. Therefore, the knowledge afforded by the present study can 
contribute greatly to the effective use of aerially disseminated pesti- 


cide sprays. 


Chapter 5S 
EVAPORATION RATE FACILITY DESCRIPTION 

5.1 General Overview 

The experimental facility utilized in the evaporation rate study 
is shown schematically in Figure 1. For purposes of illustration, the 
microscope and Reticon camera are shown above the tunnel along with the 
traverse and aerosol generator. The microscope and camera were actually 
mounted benéath thestest section jduring. testing. The tunnel consisted of 
a blower-driven closed circuit wind tunnel powered by a DC motor and 
controller. The tunnel] air temperature and relative humidity were 
altered separately with individually controlled equipment. A vibrating 
reed aerosol generator, a specially designed and fabricated long focal 
length microscope, a Reticon solid state camera and camera controller, 
and a Texas Instruments microcomputer completed the system. The evapor- 
ation rate histories of virtually any solution ranging in aerodynamic 
size from 50 - 1000 microns, ambient temperatures from 5 - 30°C, and 
relative humidities from 20 - 90% could be obtained with the use of this 
facility. The data could be accumulated with a constant test section 
velocity or with a continually adjusted velocity equal to the settling 
velocity of the observed droplet. 

The testing procedure will be described in detail later. A short 
description 1s given heretfer-clarification. First, the moving air 


stream was conditioned to the desired value for temperature and 
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relative humidity. Then a droplet of the desired aerodynamic size was 
generated and collected on a small wire (5S - 25 microns in diameter) 
held perpendicular to the airflow. When a pendant droplet was captured 
within the field of view of the microscope, Reticon camera and data 
acquisition were started manually. The data acquisition consisted of 
determining the diameter of the observed droplet at each time step. The 
following section will discuss each component of the Evaporation Rate 
Facilaty separately. 
soe Herosol, Generator 

A vibrating reed aerosol generator was developed for the purpose of 
producing droplets of a range of sizes between SO and 1000 microns. The 
vibrating reed aerosol generator was selected for use in this facility 
aS it) 15 Capable of generating targe, predictable droplet sizes and 
Decause Of ics flexipility an controlling the placement of the droplets 
generated. This system allowed the operator to place a stream of 
droplets throughout the test section. The collection of a single 
droplet on a supporting wire was facilitated greatly with the use of 
this areosol system. The vibrating reed system presented by Wolf (14) 
was used as the basis in the production of the present system. 

The entire generating system, depicted in Figure 2, was mounted on 
a traverse that allowed the generator to be moved longitudinally as well 
as laterally. The traverse consisted of four 3/4 inch diameter steel 
bars with two sets of four SPB-12 pillow blocks. A hand crank connected 
to a threaded vertical bar provided the needed vertical movement. The 
combination of these three movements allowed the direct placement of 
the generating system at any desired location within the test section. 


Four additional adjustments were contained in this generating system to 
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enable a change in the aerodynamic droplet size. Three micrometer 
adjustments modified the placement of the fluid reservoir with respect 
to the vibrating reed. The fourth micrometer changed the distance 
between the magnet causing the motion in the upper portion of the 
vibrating reed's shaft. This provided for adjustment of swing ampli- 
tude of the vibrating reed, enabling the adjustment of the generator to 
any desired aerodynamic siZe. 

The reed consisted of a 12 :mchilong stainless steel plate with 4 
thickness of 1/16 inch. The upper portion was fixed, while the lower 
end was free to move. At this lower end, a small connector was used 
to support the desired wire sizemwnicheacted as the reed) Selection of 
wire diameter was contingent upon the required droplet size. 

DPhes fluid reservoirsin -whieo tnesrecdauipratcdsconsistegsOleducloL 
Cuteine awl, X 5/40X. leainch precesot splexip lassie | WomcuCch reservoirs 
With different slot thicknesses were used in this study.) Larger droplets 
were easily generated using a reservoir with a, 1/8) inch slot, while’a 
1/16 inch slot thickness was used to obtain the smaller droplet sizes. 
The reservoir was held in place by a 3/4 X 1/8 inch aluminum bar 
approximately 13 inches long. The upper portion of this bar was connec- 
ted to the aforementioned three positioning micrometers. As these 
micrometers were adjusted, the entire bar holding the reservoir was 
moved. 

The fluid reservoir had a small diameter hole drilled from the 
downstream side through the solid plexiglass above the sldt. This 
hole terminated in the center of the slot and was used to provide a 


path between the slot and the infusion pump. The infusion pump 
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compressed a syringe, thereby generating a continuous stream of fluid 
into the fluid reservoir. A long plastic tube with the appropriate 
fittings provided the connection of the syringe and the fluid 
Resernvoyn. 

A small rectangular flat faced electromagnet, model EM-622, was 
used as the device causing the vibrating motion. The magnet was driven 
by an alternating current signal. This caused the magnet to change 
podarity,,:and thus the shaft. of -the,yreed,was,forced, to,.vibrate. ). Ai cir- 
cuit was developed to allow the input of the alternating current to the 
magnet. The alternating current was supplied by a BK Precision solid 
State model E-3108 sine/square wave generator. This wave generator 
allowed the frequency to be adjusted so as to also optimize the droplet 
generation. A General Radio 1531-AB strobe light was tied into the 
Signal generator, allowing the operator to observe the motion of the reed 
and the path of the droplets being generated. A delay was connected to 
the. strobeyso..that, both the generated droplets and the vibrating reed 
could be seen at any location of the continuous cycle. Liquid droplets 
ranging in aerodynamic size from S50 to 1000 microns could be generated 
WithetbisSeeystems 
5.3 Humidity and Temperature Control 

A wide range of temperatures and relative humidities were attained 
with the equipment installed in this facility. The tunnel air temperature 
and humidity were controlled with the aid of sensors located just 
upstream of the test section. The temperature sensor was a thermistor 
which registered relative values. via changes in internal resistance. A 


lithium chloride temperature compensated sensor was used to control and 
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detect changes in relative humidity. The exact test section temperature 
and relative humidity were monitored with a pSychometer located down- 
stream of the test section. Both the temperature and humidity sensors 
were wired into a control circuit housed in a metal cabinet. Set 
controls for the temperature and humidity were contained in the control 
circuitry and were mounted on this cabinet. The desired values of 
temperature and humidity were obtained by manually adjusting the set 
controls until the conditions were met. The controlling circuitry main- 
tained the values of temperature and humidity at the set values. If 

the temperature drifted below the required value, the circuitry acti- 
vated the heat exchangers. Similarly, if the temperature was too hot, 
the air compressor and refrigeration unit was activated and cooling 
initiated. The relative humidity was maintained in an analogous manner. 

To provide a greater control in the temperature and relative humi- 
dity, each of the systems were provided with toggle switches. The 
heating and cooling systems had two-way toggle switches that allowed 
each tobeterther -a»part of {the controfling system or compl étely shut 
off. The humidifier and dehumidifier had three-way toggle switches, the 
first two options being identical to those of the temperature controls. 
The third option allowed the system to remain on constantly, independent 
Of the sensor's. Signal’ to, the comtrollingicircuic. 

The heat used in the facility was supplied by an electrical heating 
coi that@covered the entwrerflowlarcats Thesatrewa Ss cooled usingea 
compressor, condensor and receiver system that used 45 ounces of Freon-12 
flowing through an evaporator coil design. Both systems were located 
just upstream from the test section to allow better control of the 


temperature. The control circuitry mentioned above alternately activated 
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the heating and cooling systems depending on the temperature of the wind 
tunnel air as compared with the desired testing temperature. 

The system used to maintain the humidity was somewhat diterent 
from the system used for the temperature control. When the air humidity 
was higher than the required test humidity, the air was diverted through 
a Dryomatic Mariel R-150 Dehumidifier. This dehumidifier continually 
regenerated 50% of its desicant by way of a rotating bed. The air was 
diverted from the closed loop wind tunnel by using two large, automat- 
ically controlled butterfly-type valves. One of the valves was located 
at the diverted section entrance and the other was located downstream 
in the main tunnel. Both valves rotated +90° at their center and were 
spring-loaded to remain closed unless they were forced to rotate. A 
rope and pulley were attached to the center axis of these valves and 
were interfaced to electromagnets. These electromagnets were supplied 
current by the controlling circuit, thereby forcing the valves open. 
When the humidity level exceeded the desired level, the air flow was 
diverted by closing the main tunnel valve and opening the diversion 
tunnel valve. The opening and closing of these valves was performed by 
the controlling circuitry by supplying current to one of the two electro- 
magnets interfaced with the valves. The third toggle switch configur-« 
ation mentioned earlier forced the appropriate valve to remain open 
continuously. As with the temperature system, the sensor's reading was 
used by the controlling circuitry to maintain the humidity. When the 
humidity level was too low, a water atomizer was used to raise the 
humidity. This atomizer was fitted with a heater to allow large 


amounts of water to be supplied to the air. The atomizer was located 
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in the section upstream of the test section and delivered the water 
through a four-way manifold to allow uniform humidity across the tunnel 
testaisect lone 
5.4 Velocity Monitoring 

Determination of the test section velocity was accomplished with the 
aid of a temperature compensated hot film anemometer. The sensor was 
placed in a reduced cross section area in the facility, with 10° maximum 
converging angles to prevent separation and allow accurate calculation 
of the low velocities realized in the test section. 
ovo sMicroscope 

Magnification of the image of the pendant droplet was necessary to 
simplify the acquisition of the evaporation rate data. A microscope 
developed by Dick Rossmiller Company of Boulder, Colorado was used to 
produce an image of workable size. The microscope's specifications 
included a range of magnification of 2.5X through 2500X and a working 
distance of .700 to 7.000 inches. It was mounted on a stage that 
allowed limited movement in three directions, twelve inches vertically 
and two inches both longitudinally and laterally. The back lighting for 
the microscope was provided by fiber optics. A Dolan-Jenner Fiber-Lite 
illuminator, model 170-D, was used as the light source. The light was 
also provided with movements both horizontally and laterally to allow it 
to be positioned directly over the microscope. Three leveling adjust- 
ments were also used to ensure the acquisition of uniform illumination. 
A columnating system of lenses was used with the back light to supply 
uniform parallel lighting. The inside of the microscope as well as its 
various attachments were lined with black felt or painted black to 


reduce reflection. 
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5.6 Camera and Microcomputer Description 


To determine the evaporation rates of droplets under free fall 
conditions it was necessary to acquire data on: (1) the droplet size, 
and thus allow computation of the carresponding Bereta neat ocrtty 
(2) the wind tunnel velocity, and the correction of the wind tunnel 
velocity to the calculated value. A measurement and control scheme 
of this complex operation required the use of a microcomputer. The 
Evaporation Rate Facility at Colorado State University Aerosol Science 
Laboratory used a microscope fitted with a Reticon solid state camera 
to measure the droplet size. A Texas Instruments TI-990 microcomputer 
was used to perform the necessary computations and control operations. 

93 6 1) Syst emwOperation 

A pendant droplet, held stationary in the microscope viewing field 
by a small stainless steel, tungsten or platinum wire, was imaged on a 
100 X 100 detector array within the Reticon camera. The Reticon camera 
controller provided the timing and video signal processing necessary to 
form an image on a video monitor. The resulting binary video output was 
obtained by comparing the signals from the camera to a threshold set by 
the operator. The droplet size was determined by the number of photo- 
diodes occluded above the threshold. This count was proportional to 
the average projected area of the pendant droplet. The computer read 
the data from the counters and calculated the droplet size using cali- 
bration data stored in memory and a simple mathematical algorithm. The 
settling velocity corresponding to the measured droplet size was 
calculated and compared to the actual tunnel velocity as determined by 
the air meter, analog multiplexer, and analog to digital converter (see 


Figure 1). The measured droplet size, wind tunnel velocity, and 
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elapsed time from the start of the test were displayed on the CRT 
terminal and simultaneously transmitted to the Colorado State University 
Gy DerPecomputer torastorage. 

Two modes of operation were required of the data acquisition system, 
the time measurement and the calibration mode. The typical time step 
used in the data acquisition was approximately one second. The count 
mode or calibration mode was used to determine the calibration factors 
each day. While operating in this mode, the microcomputer displayed on 
the terminal the number of photodiodes not illuminated above the 
threshold during the previous frame scan. This was done by subtracting 
the diodes illuminated above the threshold from the total number, ten 
thousand. No other calculation or control operation were performed in 
this mode. 

5.0.2 'Retacon Camera -andaControlier 

The Reticon MC520 camera was based on a solid state image detector 
formed by a 100 X 100 photodiode array (10,000 pixels). Center-to- 
center sspacing Of the pixel sewas OUfmi1CTons a heme ectroni CaciT cud try 
needed for scanning the diode array was contained within the camera. 

The RSS20 camera controller provided the camera power and the circuitry 
used in the synchronizing of clock signals for scanning. The photodiode 
scanning rate was determined from controller settings and could vary 

from 25 - 500 frames/second (250 kHz - S MHz). Camera focusing was 
achieved manually by observing the video monitor driven by the x,y,z 
display monitor drive signals from the controller. Necessary adjustments 
with the microscope were then made until the image on the video monitor 
was focused. As each pixel was scanned, the camera output to the 


controller a voltage pulse proportional to the light intensity on that 
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pixel integrated over one frame period (1/25 - 1/500 second). The 
amplitude of each voltage pulse was compared to the voltage level of 
an adjustable threshold in the controller. The threshold voltage level 
produced a bi-level digital video, the same as the pixels above the 
threshold voltage. Counting the number of true pulses, pixels above 
the threshold, for one frame scan provided a number proportional to 
Lue aropletrsice during thaterrame ‘periods 

5.6.5 Computer Configuration 

The TI-990 microcomputer components included: 
1) TM 990/101M Microcomputer Board 
2) TM 990/201 and TM 990/206 Expansion Memory Boards 
3) TM 990/302 Softward Development Board 
4) TM 990/310 48 Bit I/O Board 
S) TM 990/512 Prototype Board 

The TM 990/512 prototype board served to interface the microcomputer 
to the velocity and droplet size sensors. This board consisted of 
specialized circuits wired by Aerosol Science Laboratory personnel. 
Digital counters and their associated control circuits interfaced with 
the Reticon RS520 camera controls. An analog multiplexer and analog to 
digital converter (A/D) allowed the measurement of wind tunnel velocity 
by the computer, and a digital to analog converter allowed computer 
control of this velocity. Transfer of data from the counters and A/D 
to the microcomputer board was facilitated by the TM 990/310 48 bit I/0 
board. Software, written by Aerosol Science Laboratory personnel, 
resided in RAM on the 101M and 206 boards and in RAM and EPROM on the 201 


board. The computer operating system resided in the 302 software 
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development board. The Cyber interface and overall system operation 
was controlled by the 101M board. 

9.0.4 Calibration Procedure 

The Reticon camera was calibrated daily by determining the number 
of photodiodes occluded by a known area. The step-by-step procedure 
for calibration follows: 

1) Illumination intensity and uniformity were adjusted while 
monitoring the z axis video on an oscilloscope. | [he intensity waseset 
at approximately 90% of saturation. 

2) The data threshold in the camera controller was set at 
approximately 80% of saturation. 

3) An opaque ink stained dot of known diameter was observed by 
the microscope. The count of the number of diodes occluded by this dot 
was obtained with the computer in the calibration mode. This count was 
taken several times with the dot at random locations in the viewing 
field. The count was averaged and labeled NCAL. The incremental area 


corresponding to a single diode was determined by 


2 
_ m(Diam) 
Oo anes ANGAL (5.1) 


where Diam was the diameter of the ink dot. 

4) A stainless steel, tungsten or platinum wire of appropriate 
diameter was positioned in the viewing field. The count of diodes 
occluded by the image of the wire alone was taken several times and 


averaged. This average value is named NWIRE. 


Chapter 6 
REAGENTS 
Time history evaporation rates for 13 mixtures along with water 
were to be studied. The solutions correspond to common pesticides 
and pesticide components used today. These solutions and their 


densities an gramseper cuvic centimeters can be,seen in Table 2. 
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Table 2 


Solutions and Their Densities 


Mixture Density (2m/cec ) 
0.5 gallons of Sevin-4 oil 0.88062 


0.48 gallons of No. 2 Fuel oil 
0.02 gallons of Automate Red B dye 


4; ounce of Polyhedra Douglas Fir eos soil 
Tussock Moth Nucleopolyhedrosis 
Virus 


k gallons of Molasses 

0.7 gallons of water 

NaOH-amount sufficient to make oh 
of mixture between 6.0 and 7.2 


0.380 gallons of Sevin-4 oil 0.98770 
0.48 gallons of No. 2 Fuel oil 
0.02 sallons of Automate Red B dye 


0.98 gallons of No. 2 Fuel oil 0.86090 
0.C2 gallons of Automate Red B dye 


0.50 gallons of Thuricide 163 1.12600 
0.50 gallons of water 


4 nounds of Orthene 75S 1.02004 
1 gallon of water 
l pound of Orthene 75S 1.0428 


lL gallon of water 


2 quarts of Dow Esteron (2,4,5,-T) PeO12 24 
10 sallons of water 
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Mix Index Mixture Density (gm/cec) 
9 2 quarts of Dow Esteron 99 1.00882 


10 gallons of water 


10 4; pound of Dipel 1.05404 
4 pound of Shade 
1 gallon of water 
0.016092 pounds Rhodamine B extra 


re § 3; pound of Dipel 1.03008 
1 gallon of water 
0.016692 pounds Rhodamine B extra 


1 water 1.00000 


65 S ounces of Nalco-Trol 1.00070 
100 gallons of water 


14 10 ounces of Nalco-Trol 1.60140 
100 gallons of water 


Cha pte sa, 
EXPERIMENTAL PROCEDURE 
7.1 grestang Procedure 

The experimental test conditions used in this study are given in 
Table 1. All combinations of the three variables were obtained by 
holding two constant and performing tests for the different values of 
the third. Once the conditions of the test were selected, the wind 
tunnel was sealed and the desired temperature, humidity and initial 
velocity of the tunnel air were set and allowed to stabilize. The 
initial wind tunnel velocity corresponded to the settling velocity of 
the droplet to. bese enerated . 

Once the temperature and humidity had stabilized at the desired 
values, droplets of the given solution were generated by the vibrating 
reed aerosol generator. A single droplet was collected on the suspended 
wire within the viewing field of the microscope. Droplet diameters 
within a range of 410% of the nominal diameters were deemed acceptable 
for the tests. Data collection and wind tunnel velocity control was 
initiated with a single command to the microcomputer and continued until 
stopped by the Operator. # The researcher basedythe  decisionsto stop a2 
test by viewing data displayed on the CRT terminal. At every second 
the droplet diameter, wind tunnel velocity and elapsed time were dis- 
played on the» terminal.) sThe weneralycriteriavusedsby the researchersto 


Stop the test are: listed in Tablemsa these -eenerale test -endinomeriterica 


od 
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General Test-Ending Criteria 


The time interval of any test was to be the 


minimum of the following: 


1) The time required for the droplet to undergo 
a 90% change in original diameter. 

2) The droplet diameter becomes less than 
10 microns. 

3) The time it takes for the droplet to fall 
200m Teel. 


4) 30 minutes. 
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were contractually predetermined and represented limits that occur in 
actual spraying operations. The behavior of the individual solutions 
could also dictate separate test-end rationale for each solution. 
These individual criteria became apparent during the collection of data 
and will be discussed when appropriate. 
7.2 Adjustments in Facility Equipment 

Several adjustments of the facility equipment were required when the 
initial droplet diameter was varied. Different sizes of reed diameters, 
ranging from’25e= 250 microns, were used for the vibrating reed saerosol 
generator. The diameter of the reed depended on the droplet diameter to 
be generated and the pesticide solution. The physical properties of 
each solution varied considerably. This required the use of different 
reed) 5izes tand micrometer Settinys to generate droplets of thessame 
diameter between different solutions. In general, a reed diameter of 
approximately one half of the proposed droplet diameter was used by the 
aerosol generator. The micrometer settings were then varied until a 
droplet within the respective range was generated. 

The magnification of the microscope was also changed as the 
droplet diameters were varied. This was accomplished by changing the 
eyepiece of the microscope. During tests of the large initial droplet 
diameters, 250 and 400 microns, a 10X eyepiece and a 10X objective lens 
were used with the microscope. A 20X eyepiece was used with the 10X ob- 
jective lens in those tests which required observation of the smaller 
initial droplet diameters. This change was necessitated to eliminate 
Sizing errors that would occur by viewing the small droplets with the 


original magnifications. 
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The last required adjustment was the choice of the diameter to be 
used for the supporting wire. Three different wire materials and sizes 
were utilized throughout this study. A platinum wire with a diameter of 
17.3 microns and a stainless steel wire measuring 25 microns in diameter 
were used arbitrarily while performing the 250 and 400 micron tests. 

A 5S micron tunsten wire was used to collect the smaller diameter drop- 
lets. The reason for the change in supporting wires for the smaller 
droplet diameters was to reduce the pixel count due to the presence of 
the wire. The count of the larger wires promised to be approximately 
the same as the count of the smaller diameter droplets. This was felt 
to produce significant errors in measuring these droplet diameters. 

The use of the smaller wires eliminated this problem as well as facili- 
tating the collection of the smaller droplets. 

7.3 Time Delay 

One additional comment concerning the experimental procedure must 
be made. An increment of time existed between the actual collection of 
the droplet on the wire and the start of a given test. This time delay 
was unavoidable when using the described experimental procedure. In 
order to account for this delay, a stop watch was used to record this 
difference in time. The stop watch was started the instant the droplet 
settled on the supporting wire. The difference between the elapsed 
time recorded by the test and the time recorded by the stop watch was 
found to be the time delay. An average value of approximately six 
seconds was used for those tests in which a time delay was not recorded. 
This time delay was required in order to determine the actual initial 


d2.0 pleted tam eter, 
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1A. Transparent Solutions 


Some of the pesticide solutions used in this study were transparent. 
A problem became apparent when droplets of these solutions were observed 
with the microscope. The back lighting required by the microscope caused 
the center portion of the droplet to be too light to occlude the 
pixels.) has, glare,exceeded the threshold devell setein the Reticom 
Camera Controller, thereby producing an erroneous droplet diameter to 
be calculated. To eliminate this problem, Rhodamine B dye was added in 
small amounts to the transparent solutions. This addition increased the 
Opacity of the droplet and enabled the actual droplet diameter to be 
measured. Great care was taken to add negligible amounts of dye to 
preclude any effect on the evaporation process. Concentrations less 
than 0.2% were used. 
7.5 Data Collection and Reduction 

The values of NCAL and NWIRE, described in Section 5.6.4, were used 
by the microcomputer to determine the evaporating droplet's size at each 
time step. The pixel count due to the evaporating droplet and its sup- 
porting wire was taken at each time step when the computer was in the 
testing mode. Each of the counts was proportional to the area covered 
by the evaporating droplet and the wire. Therefore, the field-of-view 
area covered by the evaporating droplet and wire was: 

(count) AA = AREA Chat) 

The methodology used for the diameter calculation is best seen by 
referring to Figure 1 and noting that the area occluded by the droplet 
and wire in the viewing field can be divided into three sub-areas: 


1) The cross-sectional area of the droplet not common to the wire. 


pe 


2) The area common to both wire and droplet. 
aerneetotal weross sectional area of the wire. 


Also note that the width of the wire can be approximated by: 

W = (NWIRE)/100) (aa) 2/7 (7.2) 
The assumption that the wire diameter was small in comparison to the 
droplet diameter yields the following quadratic: 


2 
AREA = ———— - (W)(Diam) + NWIRE AA 
i) 


: 
(Count) AA = en 4 a Peo) WIRE. KA 


ine thesequation, count rerers to the number of diodes occluded by the 
evaporating droplet and wire. The AA and NWIRE values were determined 
Wathethe ealipracion procedure each day. Therefore, the resulting 
equation could be solved yielding a droplet diameter, Diam, at each time 
step. This diameter was used in the following equations for drag coeff- 


1eient and settling velocity to obtain the required wind tunnel velocity. 


6 


s 
Cy = 2 Pera ES, (Taken from 
Reference 12.) 
1/2 
oe 4/3(Diam) g(p_ - Pp) 
T pase (7.4) 
Dia 


The wind tunnel velocity was adjusted to this value by the microcomputer 
ar each timesstep a the value of ee or the density of the droplet, was 
assumed to remain constant throughout the test. This simplifying assump- 
tion was used to avoid the complex relation of the droplet density as a 
function of time. This relation was unknown and would require extensive 
research to determine. For this reason, the assumption of a constant 
density was believed to represent a reasonable approximation of the 


actuals physical process, 
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Chapter 8 
RESULTS AND DISCUSSIONS 

8.1 Experimental Data 

The data obtained in the study consisted of the droplet diameter 
recorded at a time step equal to one second. In addition to this dia- 
meter, the actual wind tunnel velocity was recorded. Each set of data 
was curve fitted to a third order polynomial derived to present an 
averaged but accurate result. Along with the third order equations, 
values of the average errors between the actual data and the third 
order functions were calculated. These terms, called Delta, represented 
a closeness of fit value for the third order smoothing. The third order 
functions were used over a limited range for some of the highly volatile 
solutions tested at the small droplet diameters. This limited range was 
necessitated by large errors encountered by forcing the entire test to fit 
a polynomial constant after the evaporation process had ended. For these 
Cases, the ranpe of time win ehachwene third order function was applied 
accompanies the equation. The complete set of raw data obtained in this 
study was stored on a magnetic tape. The order of tests and the required 
accessing procedures are listed in Appendix 5S. 
o-c blimination-or fests 

8.2.1 Solutions of Low Volatility 

Among all of the 14 pesticides and pesticide components, two dist- 


inct groups of solutions were prevalent. The first group contained 
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water as part of the total mixture and had a tendency to evaporate 
relativelysfastwasl hewsecond! group} esolutions ly grand) 4), ycontainedsoil 
and evaporated at a much slower rate. Close examination of the data 
obtained for the second group of solutions revealed two notable tend- 
encies. First, when working with the larger diameter droplets, 250 and 


400 microns, the total change in droplet diameter was less than 5 


of 


The 
initial tests were performed at the extreme conditions felt to produce 
the greatest amount of mass transfer--i.e., high temperatures and low 
humidities. Since the change in diameter at these critical conditions 
was small, the change in diameter using low temperatures or high humi- 
dity values could also be assumed to be small. Therefore, tests 

that were to be performed at low temperatures or high humidities 

were’ eliminated. 

Examination of the smaller diameter tests, 50 and 100 micron initial 
diameters, revealed the second trends. Tests were performed using these 
solutions over wthe time: intervals required for the droplet to fall 300 
feet. Examination and comparison of the tests performed at 20 and 60% 
relative humidities showed little or no difference. It was concluded 
that raising the relative humidity to 90% would produce results not 
achieving the minimum criteria of Table 3. Therefore, the tests were 
eliminated with 90% relative humidity at all temperatures. Tables 
showing the rationale used to eliminate the tests for solutions l, 3, 
and 4 can be seen in Appendix l. 

6. 2c ooLut 10s: Off Haph *Volatil ity 

Time rate evaporation histories were to be obtained for solutions 
13 and 14 which contained different concentrations of Nalco-Trol and 


water. After acquiring data under a variety of test conditions, the 
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Nalco-Trol additive was found to have no significant effect on the 

rates of evaporation of the droplets. The result was that the time rate 
histories were essentially identical to the previous data determined for 
water. To ensure this result, a new solution was mixed which increased 
the amount of Nalco-Trol added to the water by a factor of ten. Examin- 
ation of the data corresponding to this new solution produced similar 
results--the Nalco-Trol did not detectably retard the evaporation 
process: 

Another test was developed and performed to support the results 
obtained using the evaporation rate facility. Approximately 10 cubic 
centimeters each of water, solution 13 and solution 14 were measured 
into 100 ml beakers. These beakers were set side by side and weighed 
at prescribed time intervals. The mass lost by free evaporation from 
the beakers was then compared. The mass loss of each solution was 
normalized to a percentage by dividing by the total) mass. “These percent 
mass losses were then plotted as a function of time and compared (see 
Appendix 2). The comparison showed a negligible difference in the loss 
of mass between the solutions at all times. Therefore, the conclusion 
was sustained--solutions 13 and 14 evaporated at the same rate as water. 
The remaining tests were then eliminated for solutions 13 and 14. 

It should be mentioned that the addition of the Nalco-Trol did 
Significantly affect the generation a kane droplets. The Nalco-Trol 
additive greatly increased the surface tension or "elasticity" of the 
solution. The result was that the solution would stretch considerably 
before the filament of fluid would relax and form a droplet. This indi- 


cates that the Nalco-Trol additive may serve to reduce the number of 
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fine droplets formed, thus reducing drift loss and perhaps narrow the 
resulting size distribution of droplets formed by presently used atomi- 
zation systems. 
8.3 Small Droplets of Volatile Solutions 

Cx oT ee Lene tn or Lests 

Tests performed on solutions 2 and 5 - 11 that contained large 
amounts of water were seen to evaporate rapidly using the small initial 
droplet diameters. Most of the change in droplet diameters occured in 
thestime Tequiredetor crystallization to begin. | In tests involving 
these solutions, a relitively constant diameter particle remained on 
Die wi resallereathe shirstsminute, of thestest.. To insure no additional 
change would be experienced, each test was performed over a period of 
ZOUSSeCcOnds seinmallecases none 01 the; test-ending criteria of Table 3 
were met during this interval. However, it was assumed that the 
evaporation process had been completed and the tests could be ended. 

82522 Approximation of Tests 

An empirically calculated technique was used for solutions 2 and 
5 - 11 with initial diameters of 50 microns and additionally an initial 
droplet diameter of 100 microns for solution 5. The rapid evaporation 
and difficulty found in generating and collecting small droplets required 
for these tests necessitated the application of this technique. The 
parameters, assumptions, and calculations used in this approximating 
technique can be found in Appendix 3. 

The technique used to approximate the small diameter tests relied 
on four assumptions. Application of this technique also required the 
evaporation process of a given solution to be divided into three separate 


stages: normal.evaporation, transition and crystallization. The first 


40 


stage, normal evaporation, was the portion of the process in which most 
of the mass was transferred. It was this stage of evaporation that was 
determined using the first assumption. This assumption stated that a 
water-based solution evaporated at the same rate as water until the inert 
material contained in the solution began to crystallize. Direct observ- 
ation of the present data supported this assumption. This tendency 
remained consistent for all temperatures and relative humidities. 
Several authors, including Ranz and Marshall and recently Picot, 
Chitrangad, and Henderson (9) have reported data that also supports this 
assumption. Therefore, the first stage of evaporation can be approxi- 
mated vy using the water data obtained by this study. The second-by- 
second droplet diameters used in the approximated tests were taken from 
the water data with the appropriate environmental conditions. This was 
continued until the crystallization stage began. 

Determination of the droplet diameter at which the crystallization 
started (called the critical droplet diameter) required the application 
of the second assumption. This assumption stated that the crystalliz- 
ation droplet diameter depended on the solution concentration and initial 
droplet ‘svze-alone. The ratio of the volumes of the critical’ droplet to 
the initial droplet, as stated in assumption three, should then be a 
constant value for any initial size. THis was to remain consistent for 
any temperature and relative humidity combination. The data obtained 
for the initial droplet diameters of 100 microns were used to obtain an 
average ratio of volumes for critical to initial droplet diameters. 
Therefore, the critical diameter for an original diameter of 50 microns 


could be approximated. This critical diameter used in combination with 
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the first assumption and the appropriate water data was used to approx- 
imate the first stage of evaporation. 

At the point of crystallization the rate of evaporation reduced 
Signatacantly sey hism@ratesreductiongwas catresultqof the decreasing 
amount of volatile material present. As the evaporation progressed, 
the volume of the volatile material was continuously reduced and more 
inert material crystallized. As more and more material crystallized, 
the rate of evaporation reduced accordingly. The process of reducing 
the rate of evaporation gradually was seen to be the cause of a curved 
transition between the first and second, or crystallization, stages 
of evaporation. This curved transition represents the transition 
stage of evaporation and will be explored later. 

In approximating the second stage of evaporation, a fourth 
assumption was needed. This assumption was simply a direct observation 
obtained by- examining several tests. The second stage of evaporation, 
the process in which the critical diameter reduces to a final constant 
diameter, was seen to be closely approximated by a constant slope line. 
This slope was seen to be relatively invariant to changes of temperature, 
humidity and initial droplet size. Therefore, once a value for the 
slope had been found, it could be used for all droplet sizes and environ- 
mental conditions. An average value of this second stage slope was taken 
directly from the 100 micron tests for each solution. The second stage 
slope was continued until the evaporation process had ended, or the final . 
droplet diameter was obtained. The final droplet diameter, like the 
critical droplet diameter, was considered to be dependent upon the 


original diameter and solution concentration alone. A ratio of the 
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final diameter to original diameter for each solution was found by 
averaging the values given by the 100 micron tests. A volumetric con- 
centration relationship could also be used, however the volume of inert 
material must be known to implement this approach. 

Combining the approximation techniques for the first and crystal- 
lization stages produced an approximation for the entire evaporation 
process. The application of these two techniques required the use of 
the data obtained for water and the given solution with initial diameters 
of 100 microns. Once the parameters used by this technique were 
obtained, two sets of diameter-versus-time approximated data were deter- 
mined and combined. The transition stage of the curve was then found 
by applying a numerical smoothing subroutine. The subroutine changed 
the sharp transition between the first and second stages to a anoeEh 
gradual curve. This was consistent with the gradual transition between 
the twoustagesvas Seen in actuals tésts: 

Appendix 3 also contains tables listing tests and calculations 
required to determine the various parameters used in the technique. 
Statistical values for the mean absolute and probable errors obtained 
inethisete chniquésaretgivenstor téach> solutions ss Plotsrofrathe critical 
versus original diameters are also shown. Summarizing the errors pro- 
duced by the approximation technique reveals all are less than 9%. The 
maximum change in the critical diameter corresponding to this error 
would be = 2 microns. This is equivalent to a volume difference of less 
than 25%. Therefore, application of this technique enabled data that 
could not be obtained experimentally to be calculated empirically with 


conf idéncemin®theraccuracysofepredictions. 
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5. ovo, 00 Microm ests for Solution 12 

Included in the tests using an initial diameter of 100 microns were 
the data with an or*zinal diameter sf 50 microns. Since solutions 12, 
water, was a homogeneous solution, the evaporation should continue uni- 
formly throughout the test until the substance has completely evaporated. 
The dye needed to view and perform tests for water at 500 and 250 microns 
was not required for the small diameters. Hence, the tests performed 
With the Original diameter” of YO0*microns *representeéd’ pure water:* “The 
required data with an original diameter of 50 microns could then be 
taken directly from the 100 micron” tests’ performed. 

8.4 Simuzation’ or Terminal’ Velocities 

The experimental facility used in this study was to determine the 
Tales iess eransrer or drop. etsital lino at therr *termuial wélocities . 
Droplet diameters were measured at each second and their corresponding 
settling velocities calculated with equation 7.4. The wind tunnel 
velocity was then adjusted to the terminal velocities by the microcom- 
DUGere once ae err ie se Cts >. or. Tew actual Test velocity’ was 
recorded along with the droplet diameter and elapsed time. 

Figure 3 shows a comparison of the calculated and the actual test 
VElOCIt1es 101. 2 representative test... Examination of this’ plot#reveals 
that’ thei two velocities were very close at ‘each step.’"\There was a 
Slipmiecey tation ieenee Ve lLOGLULesratla het Deginning of tne test -© This 
Waseem rect. tesult Or pie time required for the “computer “to initiate 
Coltrol or tie velocity, “Ihe test velocity*also-tended to be a little 
larger than the calculated velocity. This indicated that there was a 
time delay between the determination of the droplet diameter and the 


adjustment of the wind tunnel velocity. ‘This resulted in’ ‘the test 
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wind velocity lagging behind the calculated velocity by approximately 
one second. However, it is assumed that the differences between the 
actual test and calculated velccitiesS were negligible. Therefore, the 
terminal velocities of the droplets simulated in this facility were 
Pelt itosbe accurate: 
8.5 Evaluation of Data 

Sgr aly Water Data 

Of all the solutions examined in this study only one, solution 12, 
can De explored in the traditional manner presented by Froessling, Ranz 
and Marshall, and others. The determination of an equation which 
presents the Sherwood number as a function of Re and Sc required the use 
of several little known physical properties. Values of the diffusion 
coerficient and the saturation vapor pressure are not easily obtained 
POL alt SOLUTIONS seal iese physical properties have been explored at great 
extent for water, solution 12. Therefore, this solution was used to 
determine a relationship between the Sherwood number and Re and Sc. This 
allowed direct comparison of the resulting relationship with those pre- 
sented by previous authors. 

The parameters and equations required to obtain an equation relating 
Sh with Re and Sc are given in detail in Chapters 2 and 3. Equation 2.6 
must be used to determine values of the Sherwood number at every time 
Stop we aiica t0bene vandwoc,, civen by equation 5.1, must also be deter- 
mined at each time step. The physical properties required to determine 
these values are shown in Table 4. 

ifmOrac Veto uuseecdUatione..0 tne temperature at the droplet surface 
must be known. Ranz and Marshall, Kinzer and Gunn, and Fuchs all state 


that this temperature is basically equal to the ventialted wet bulb 
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Table 4 


Physical Properties for a Water 


Droplet Evaporating into Air 


Dry Saturation Vapor 

Aire temperature.) DiffustonsCoertfi cient gAir Density Pressure 
(se) tent eee) (ae) (bars ) 

-5 02210 0.001088 0.0041 

0 07246 0.001068 0.0061 

5 O7222 0.001049 0.0114 

10 Oa227 0.001031 Oe0LZ5 

ES 0.234 0.001013 OSOL7 1 

20 0.240 0.000995 0.0234 

50 0.254 0.000962 0.0250 


Diffusion coefficients taken from Davies. 


Values for the dry air density were calculated 


from the ideal gas equation. 


Values for the saturation vapor pressure were 
taken from the Steam Tables of Van Wylen 
and Sonntag (11). 
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temperature. Experiments have shown that this is true within about 
fe) 
O22 J* Thus, thessurface temperature of the droplet was taken as the 


wetibubo ?t emperatiresvelhe*vapor ‘pressure ‘at’ the’surface, “P_, was taken 


0” 
to be the saturated vapor pressure at the surface temperature. The 
free stream vapor pressure, P|, was determined by multiplying the air 
humidity times the saturation vapor pressure at the air temperature. A 
geometric mean value for the diffusion coefficient was used in this 


equation, as suggested by Fuchs. This value was obtained using 


equation 8.1. 


DPSs 7D > =) (8.1) 


The denominator of equation 2.6 could then be calculated. At each 

Step the only change in this value would be the radius of the droplet. 
The numerator of equation 2.6 was determined with the data obtained 

in this study. A higher order central difference method was employed to 

find the time rate of change of the diameter at each step. The data was 

first smoothed by a spline and then equation 8.2 was used to determine (2) 

c 

This mass derivative was calculated at each time step and divided by 

the denominator to obtain a value for the Sherwood number for each second. 
Most authors have related the Sherwood number as a function of Re 

and Sc. The most commonly published of these relationships has been 


pet “gah! te The value of get de 


1/2 


ture uSing water. Re 


remained constant for a given temper- 
varied as the droplet diameter decreased. 
Because the described facility simulated terminal velocities, both a and 
Ve varied in the equation for Re. This value was calculated at each 


SeconGg sana sed to wptain a value of Gua red! “sc l/3. 
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Values for Sh and C were determined for solution 12 using tests with 
every variation of temperature, relative humidity and initial diameter. 
Examination of tests using the high humidities indicated that at 20°C 
and 15°C the actual humidity attained was somewhat less than the desired 
90%. The reason for this discrepancy was the difficulty in sustaining 
the required amount of moisture in the air at these higher temperatures. 
Great difficulty was found while attempting to stabilize the humidity 


at high temperatures. The stabilization of these high humidities was 


also hindered with the small diameter droplets. The wind tunnel 
was operated at low velocities in order to match the terminal velo- 


cities of these small droplets. These low ventilationjvelocities 


resulted in slight alterations in the wind tunnel air humidity. Because 
the velocities were so low, a consistent value of relative humidity was 
not easily maintained using the high humidities, 90%. For these reasons, 
none of the tests requiring 90% relative humidity were included in the 
analysis tO =find olieas aeLunction, or Cc. 

The values of Sh and C were plotted over the ranges obtained in 
this study. A linear regression was performed over these ranges of Sh 
and C. The plot and resulting regression equation are shown in Figure 4. 
The equation representing the results obtained by this study is given 


in equations. 


a 
Sh tls. 750.4 02555 Reese. (8.3) 


A significant amount of scatter was shown by the plot and the 
value. given. for. the, average;percent epror. sihis ’seatter was a result of 
the variety ot initial conditions used. Four temperatures, two humidi- 


ties, and three initial droplet diameters were used to obtain the 
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results. The regression equation given correlated to the data fairly 
well. A correlation coefficient of about 0.8 was found and indicated 
a relatively close fit. The outer lying data points shown in Figure 4 
were produced by actual test conditions different from those recorded. 
The outer points lying below the regression line correspond to low 
Ae O 

temperature and low humidity tests, or tests performed at S C and 

20% humidity. Similarily, the outer most points above the line were 


obtained with tests using the higher temperatures and humidities. 


The representative equation, 8.3, agreed well with previous 
results. The mass transfer coefficient, §8=0.535, was very close 
to those obtained in other works. However, the intercept term, a, 
obtained in this study differed from those reported previously. 
Previous authors have determined this intercept value to be equal to 
the theoretical lower limit of Sh=2. Ranz and Marshall determined 
an intercept of 1.79 experimentally using water. The authors believed 
that the low value was caused by neglecting the change in vapor 
concentration at the surface of the drop. Beard and Pruppacher have 
also reported an intercept below the theoretical lower limit. However, 
the authors described the convective mass transfer with two equations 
for two seperate ranges of Genes! “sca A straight line was drawn 
from the last data point to the lower limit of Sh=2. In this way, the 
data was represented in two ranges in order to conform to the 
theoretical requirements as well as the data obtained experimentally. 


the presentestudy estabilished@an vinterceptvofeasliy 55.) elhis 


value is below the theoretical lower limit of Sh=2, shown in Figure 4 


by the dashed line. Limited data were obtained in the lower range of 


aye 


C. Therefore, an actual description of the regression curve in this 
lower region was not possible. However, the trend of experimental 
data justified the continuation of the regression line to C=0. 
Values below the theoretical lower limit of Sherwood number equal to 
two could be obtained using the model representing the present 
experimental data. This indicates that a discrepancy exists between 
the theory discribed in Chapter 2 and the experimental data aquired 


gn thasestudy. 


Cog tects ule: COfcilanges il Experimental Conditions 

Wafewcitects Or relative changes in temperature, relative humidity, 
@ioiiitia Gropp let aiianeler cal easivy be Seen with the use. of pilots. 
Plotting aces combination of condition and comparing them would require 
a ¢rear amounc+on time ald space. FOr thiS reason, *a representative set 
of conditions were chosen. Plots comparing the relative effects were 
performed by holding two of the conditions constant and varying the 
third. “Therefore, three plots of multiple curves were obtained for each 
solution. The elimination of some of the tests reduced the amount of 
curves on some of the plots. The representative conditions were set at 
208G? 20% relative humidity, and an initial droplet diameter of 100 
microns. Again, this condition represented the highest. rate of 


volatility. 


we 


In order to allow the direct application of the data obtained, the 
cumulative percent of mass lost was plotted as a function of time. In 
order to determine the actual percent loss of mass at each time step, 
the original droplet diameter needed to be determined. This was found 
by using a relationship of the droplet's surface area at constant vel- 
ocities: At™a constant velocity, the slope otathe surface sareagyersus 
time curve can be seen to be constant. 


Ee -2 


es : 4 
ae 5 (Cc, Ge) constant (8.4) 


Pp 


During the ‘delay times betweenethe collection of thesdroplet and the 
Start of the data acquisition, the wind tunnel velocity remained constant. 
therefore equation = $74 = couldsberused to rindache originaledroplet 
diameter. The constant slope of the surface area versus time was approx- 
imated from the first three data points. This slope was then multiplied 
by the delay time and added to the surface area of the first data point. 
The resultant initial surface area was used to find the initial droplet 
diameter. The percent loss of mass at each time step can then be cal- 
culated using the initial diameter with the diameter at each step. 

Again, the assumption of a constant droplet solution density must be 
made. As stated previously, this assumption is incorrect, but allows a 
close approximation. Using the assumption of a constant density of the 
droplet, the cumulative loss of mass at each time step was found by 


equation 3.5. 


percent mass loss = ea - cane (8.5) 
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where egies is thevinvtialedroplet-radius. Solving equation 8.5 for each 
time step allowed the curves of mass loss over time to be obtained. 

As mentioned in section 8.2, two <istinct groups of solutions can 
be seen -- cence. and oil-based pesticides. Generally, the sol- 
utions in each group behave similarly. The results of the procedure 
given above will be examined in these general groups. Distinct indi- 
vidual characteristics of each solution will be detailed during these 
discussions. 

The plots of percent loss of mass as a function of time given in 
this study can easily be used for direct application. Knowing the 
temperature and relative humidity of the surrounding air and the average 
droplet size produced by a spray nozzle, allows the operator to determine 
the approximate total loss of mass occurred over a specific time. The 


percent 0ss Of Mass for a Single droplet as given by this study: could 


©) 


be used to approximate the percent of mass lost in a spraying operation. 
ThessOray ine werent of the operation and the terminal settling velocities 
DE tite arOp lets can pesusea to determine the time antérval of the oper- 
ation. The average droplet diameter produced by a spray nozzle could be 
used with the appropriate plot to give an approximate percent of mass 
lost for the given time interval. This could then be used with the 

total starting mass to determine the amount of pesticide remaining. Thus, 
the mass ‘of pesticide’ reaching the intended’ target is approximated. This 
result would be conservative as the effects of aircraft wake and atmos- 
pheric turbulence have not been considered. A more involved procedure 
could be used to find a closer approximation if the droplet distribution 


produced by the nozzle was known. Summing the loss of mass from the 
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different droplet sizes would yield astotal ‘percent, or mass@iose. 
Again, the pesticide remaining could be determined. This remaining 
mass of pesticide would give an idea of the effect on the insect pop- 
ulation produced by the spraying operation. 

Using combinations of the resulting plots of this study would ai+ow 
the acquisition of the percent of mass lost for conditions Tiot shown: 
The relative differences between curves with varying temperatures, numi- 
dities, and) droplet srzes can be, used «toredetermine the percent, 10ss jor 
mass of a given condition. A simple example will illustrate this point 
clearly. The percent of mass losteafters five: secondssotssoOlUurloNe1o. 
20°C and 60% relative humidity with an initial diameter of 240 microns 
15 desired. “Using Figure 225" the percent 10ss Ofemasse0s wine sane rome 
with a 20% relative humidity was about 20% after five seconds. A drop- 
let whose diameter was 100 microns in the 20°C and 20% humidity air Tost 
about 75% of its mass after five seconds. From Figure 44, the same 
droplet with a 60% humidity lost only 40% of its mass after five seconds. 
Using these three conditions, the mass lossvot, tne s240Pmicron droplet 
at 20°C and 60% relative humidity after five seconds can be approxi- 
mated by: 


Dv 
mass loss (240 um, 20°C, 60%) = ae (0.40) = .106 = 10.6% 


An actual test of a droplet with a diameter of 240 microns at 20°C and 
60% relative humidity lost 10% of its mass after five seconds. There- 
fore, this procedure will allow a close estimate of the percentage of 


mass lost for conditions not included in the plots. 
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8.5.5) oO UNEIONS of. Low Volatility 


Solutions 1, 3 and 4 were included in this group of solutions. 
Each was composed of varying amounts of Sevin-4 and No. 2 fuel oil. 
The relative effects caused by changes in droplet diameter, temperature, 


and relative humidity can be seen in Figures 5S - 13. 


rc 


8.5.5.1 Change in Initial Droplet Diameter 


Figures 5 - 7 show the mass loss of solutions 1, 3 and 4 respective- 
ly for four different initial droplet diameters. The actual initial 
diameters are also given. The rate of mass loss can be seen to increase 
as the initial diameter decreases. At successively smaller ahs: 
diameters the rate of percent loss of mass increases. Figure 6, corr- 
esponding to solution 3, shows a variation to this observation. The 
two larger droplet diameters have near identical curves for very differ- 
ent initial diameters, 420 and 240 microns. This may be a result of the 
increased amount of Sevin-4 oil used in this solution. Figure 7 also 
shows near identical curves for the two larger droplet diameters. How- 
ever, thesé) initial diameters are fairly close, 330 and 280 microns, so 
that) néarevaenticalvcirves could be expected. In thesée™figures, the 
curves corresponding to an initial diameter of 50 microns have a great 
deal of scatter. Small variations in the droplet diameter measurements 
become more evident using these small initial diameters. Therefore, the 
resulting curves of percent loss of mass versus time show some scatter. 
A smooth curve was drawn through the points to represent the data. With 
the solutions of low volatility a decrease in the initial diameter from 


about 400 to 50 microns will increase the percent loss of mass from 


approximately four to close to 20% after 50 seconds. This is an increase 
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Of agiactoreos diveyansthe,percent.mass.loss for the, range of droplets 


used an) thi srstudyt 
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3a Sai 2ee Change, in, Air, Temperature 


The effect due to the change in air temperature for solutions 1, 3 
and 4 are shown in Figures 8, 9 and 10. As expected, the rate of mass 
loss increases as the temperature of the medium increases. In these 
figures, the initial droplet diameter accompanied the air temperature. 
In some instances this may become a factor in examining the plots. In 
each of the plots, the curves of increasing temperature generally 
increase gradually. The percent of mass lost using 5°C was about 3% 
after 50 seconds. For the same length of time, the 20°C curve loses 
approximately 9%. The percent loss of mass increased by a factor of 
threes forvthe temperatures used.in,this, study. 

Seo teint Change, in Alt, Relative Humidity 

Provumesm ble 1) display,the effects of change in, relative, humidity 
in the mass loss of droplets. of solutions 1, 3 and 4. Again, the dia- 
meter of the initial droplet was also given. These plots show that 
there are negligible effects on these solutions due to a change in 
humidity. This was given earlier as a criteria for eliminating tests. 
In each figure, the differences between the curves of varying humidity 
were small. An explanation for this can be shown by examining the 
equation for stationary evaporation, equation 2.4. The driving force 
in this equation was the difference of the vapor concentration between 
the droplet and surrounding air, equation 8.6. 


M P 


Ps 
DF. =—= (5 - Je CC, - C) (8.6) 
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Quantities for the vapor pressures, or concentrations, are not known 
for these solutions. However, since the vapor of the evaporating drop- 
let was basically composed of oil, the addition of water in the form of 
humidity should have little effect on the concentration gradient of this 
vapor. Therefore, the entire evaporation process using droplets of these 
solutions should change very little for changes in air humidity. This 
relatively small change applies to the limited ranges used in this 
study. Over longer periods of time; the effective differences may 
become significant. 

S064 DO LULLONS Of StL one Vota Gal 2 Gy 

Solutions that contain water as the basic component are included 
in this group.” Solutions Z-and#o= 14Pbelongsin thistnigh volatimacty 
group. Figures 14 - 46 show the effects on the percent mass loss as 
a function of time due’ to changes in’ the experimental conditions...) It 
should be noted that each of the curves was continued past the start of 
crystallization. This was done to show the different stages of evapor- 
ation described earlier. Examination of all the plots revealed that 
the critical point, detailed’ in “S.35'2) “wasian fact ‘dependentuomuite 
initial diameter and solution concentration alone. In nearly all the 
plots, crystallization Starts “at*the*same evel of percentyloss af mass 
for the varied initial conditions. This was a basic assumption used in 
the approximation of the 50 micron data for these solutions. 

8.5.4.1 Change in Initial Droplet Diameter 

Figures 14°= 24"display thereftectsvor chang imps thesanitaaladroplet 
diameter for each of the solutions. As with the low volatility solutions, 
the smaller droplet diameters lose mass atgjawfaster rate than thes larger. 


The curves used for the 50 micron initial diameter were obtained from the 
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approximated data detailed in section 8.3.2. Plots describing solutions 
13 and 14 contain only two curves. The elimination of the other tests, 
as mentioned in section 8.2.2, reduced the number of curves. Each 
figure shows the individual solution's dependence on droplet size. How- 
ever, all the plots evaporate at a rate similar to water. This was a 
basic assumption for the approximation of data in section 8.3.2. 
809.4.2— Change vin Air Temperature 

The effect on the percent loss of mass due to changing the air 
temperature for the volatile solutions can be seen in Figures 25 - 55. 
Examination of the different curves, corresponding to the different air 
temperatures, revealed the higher temperatures produce faster rates of 
evaporation. The effect of different initial droplet diameters will 
also alter these curves, therefore these diameters were given for each 
curve. An even transition from curve to curve can be seen in most of 
phe plots Glrregularities are shown inesome or thes igures es inese 
could have been caused by slight errors in the temperature and relative 
humidity readings for the=test. “For example, srigures«26me2/. 29, ol 
and 33 describing solutions 5, 6, 8, 10 and 12 show a larger difference 
between the curves for 5°C and 10°C than between any of the others. 
This indicates that the actual humidity may have been higher than 20%. 
Also, Figures 28 and 30 describing solutions 7 and 9 both show near 
identical curves for the 5°C and 10°C data and the 15°C and 20°C data. 
This may be attributable to an error in the data. Interpolation between 
the 5°c and the JORG data should give an accurate approximation. Gener- 
ally, the trend caused by the change in air temperature change can be 


Seen) with these plots. 
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8.5.4.3 Change in Air Relative Humidity 


Figures 36 - 46 show plots of the percent mass loss as a function 
of time for each solution over a range of air relative humidities. Each 
solution was described by the different figures. As the humidity 
increased, the rate of mass loss was seen to decrease. This effect gen- 
erally seemed to get more pronounced as larger amounts of moisture were 
added to the air. This observation was very much different from what 
was seen with the low volatility solutions. The reason for this differ- 
ence in the high volatility solutions can be seen by recalling the 
driving force of evaporation, equation 8.6. The vapor concentration at 
the surface of these droplets should be similar to the vapor concentra- 
tion of water. This is because the basic component of these solutions 
was water. AS more and more water is added to the environmental air, 
the C. term should increase. The total effect would be to decrease the 
cermepiven vy so.0 anid, therefore, the entire evaporation. 

8.5.5 Summary of Mass Loss Plots 

A general summary of the results given in sections 8.5.3 and 8.5.4 
will give an overall view of the effects of changing each initial con- 
dition. The individual effects on each solution can be determined by 
examining the respective figures. The physical process will be 
examined so as to determine the cause of the different effects on mass 
transrer. 

Rearranging the equation for the Sherwood number, equation 2.6, 
allows a value for the rate of convection mass transfer: 

($) 5 (2) (EEN Cr.) 
2 ae 


c : Tp 
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Assuming a functional form of Sh can be determined for all solutions such 
thateShi= ShiiRe Sc) 35 s7ebecomes: 


see ee =) (8.8) 


R et 


(S@) = Sh(Re,Sc) { 
es 0 oo 

Equation 8.8 defines the theoretical rate of transfer of mass from 

droplets. The relative effects on this value due to the changes of 

conditions can be examined with this expression. 

Of all the changes of initial conditions, the change of initial 
droplet diameter had the greatest effect on the mass transfer. As the 
initial diameters decreased (an increase in the surface area to volume 
ratio), a larger percentage of mass was lost for identical time inter- 
vals. The cause of this effect can be seen by examining the initial 
values of the masses. Because the initial mass of the smaller diameter 
droplets was itself small, small changes in the diameter as the droplet 
evaporates corresponded to large percentages of its mass. The (2) 
value from equation 8.8 can be shown to decrease for smaller eee: 
droplet diameters. This effect was doubled when using freely falling 
droplets because Re decreased as the diameter decreases. This indicates 
a reduction in the rates of mass transfer. However, this reduction is 
caused by the continuing decrease in the amount of available droplet 
mass. 

Changing the air relative humidity caused the next greatest effect 
on the mass transfer in the high volatility solutions. As discussed in 
section 8.5.3.3, the addition of moisture into the air had negligible 
effects on the evaporation of the oil-based solutions. The effect of 


increased humidity on the water-based solutions was discussed in detail 


in) section) 8.5.4.5." “Since the driving forcesterm,esiven by eougtaones..o, 
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was seen to decrease as the humidity increased, the rate of mass trans- 
fer, equation 8.8, of these solutions'must also decrease. Using identi- 
Gal instealedropletsyinadifferenteaairrhumidities; thempercentrioss*of 
masS can easily be Seen to be larger for the, droplet in'the»lower 
humidity. An example using 20% and 60% relative humidities is shown in 
Figure 47. 

The last change of testing conditions was the varying of the air 
temperature. The effects of this change were relatively small compared 
to those Gausedypyechanging therinitialijdroplet  size.ae,Changingsthe 
air temperature affected several terms in equation 8.8. Values of 
Résuses Ds Pop Re. and Ty weresalter edabyeherchangeciny ainptempeyature, 
T- The values of Re, Sc, and D should vary slightly with the changes 


in air temperatures used in this study. Realizing the temperatures, 


T) and T_, were expressed in units of degrees Kelvin, the range from 


O Lape Gena ee 
5°C to 20°C becomes insagniticantyalisayenlne values of the vapor 


pressures should vary little more than D with temperature, see Table 4 
for an example, Therefore #*changes 1ntalri temperatures have only slight 
effects on the mass transfer as compared with the other changes. The 
exception to this was the effect of varying humidity on the oil-based 
solutions. However slight, the rates of mass transfer were decreased as 
the temperatures were lowered. As seen in Figure 47 with different 
humidities, it can be shown that the percent loss of mass also decreased 


as the temperature decreased. 
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Figure 47 


Calculations Relating Rate of Mass 


to Percent of Mass Loss 


Given that- rate of mass ctranster=is 


greater tor lower humidity values 


fe ; 4 (* =) 
NG ETE aie 60% 


The time intervals are the same 


20% 
Dividing by the initial mass 
m m 
sea A ole GIST: 
Theserare percent Joss "0r Mass values 


(percent loss of mass) (percent tossear MASS) 9 
G2 


20% 2 
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S.0 Estimation of Error 

e.Gole Error pue to Supporting Wire 

Experiments using pendant droplets have been examined with extreme 
Caution by previous authors. The main concern with such experiments has 
been theyerror introduced: by the wire mised thotisupport the droplet. The 
heat, SUppized to the droplet by the wire has been the subject of most of 
Opis-Graticism: This section will attemptiito show that the error intro- 
duced by the supporting wire in this study was negligible. 

The heat transfer of a ventilated droplet is an analogous procedure 
Soetneemassetranster of tne same system. The effects of forced convec- 
tion on this heat transfer are contained in the Nusselt number, Nu. 
Slmcemionesomysicaliprocesses Of leat and mass transfer aré simijar, the 
Nusselt number and Sherwood number should have similar empirical ex- 
pressions. Therefore, an expression fon the Nusselt) number as a function 
of Re and the Prandtl number can be given. 


V2 ys 


Niele 5/5. Oa. Re Pr (8.9) 


Pod? tones. canbe used to determine values of Nu for any droplet size. 
The Nusselt number was defined as a wind factor in equation 2.7. 
=) 
R on 
Ce) : 


Si yp eee) Po boar) 
toes ~ Qnak(T-T,) 


at itdm \e 
(Ge), 


dm ; 
This expression can be rearranged and a value of (:22 | determined for 


jas 
= 


| 


Qu 
ct 





ee 


€ 
anyearop letmtaiiing treely through air. The ey térm represents the 
c 


rate of heat transfer of a ventilated droplet. Comparing the heat con- 
ducted through the supporting wire to this value will reveal the rela- 


tive amount of error introduced. 


106 


The heat conducted through the supporting wire can be calculated 
easily using the general equation for the heat flux: 


Cle gels } 
dQ Mee Roe SOR 8 ea (8.11) 


dt dL L 

where L is one half of the length of the wire. As described earlier, 
three different types and sizes of wires were used. This analysis will 
include the 17.8 micron platinum and the 5 micron tungsten wires as they 
were used most frequently. Applying equation 8.11 with the most extreme 
of the environmental conditions, 20°C and 20% relative humidity, values 
of the heat flux can be calculated. These values are compared with the 
values obtained from equation 8.10 in Table 5. It should be noted that 
the smaller wire was used when testing the smaller sized droplets. 
Thererore, the heatatlux. values ca lculateds forstheuwaresawere compared 
to the heat transfer values of the droplets with which they were commonly 
used. The largest percentage obtained by this analysis was 3.4%, which 
occurred at 50 microns using the small wire. It is believed that such 
a percentage 1sinot Signiticant. “Ihissadditioneot heats woul decausesda 
rise of the droplet temperature of approximately 1%. This is about 
the same as the error expected in measuring the air temperature. There- 
fore, the heat conducted by the wire is negligible as compared with 
the heat transfer of the pendant droplet. 

8.6.2 Estimation of Experimental Error 

An error analysis was performed for this study to determine an 
expected error of the results. The step-by-step analysis can be seen 
in Appendix 4. This analysis was used to determine expected errors for 
Sh and C for the ranges used in obtaining the plot shown in Figure 


4. This, in turn, should give a good indication of the actual 


Droplet 
Diameter 
(microns) 


400 


200 


100 


50 
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TaolLe.s 


Heat Conduction Through Supporting Wire 


Nol = io WSS 2 OeSS5 aes? 


Pr = 0.72 at 20°C from Davies 


am oy eet 
Nu es dt dt/im . Percentage 
(cal/sec) (eal/sec) 
-4 -6 
Amo? 3.80X10 3 S650 0.010 1.0 
J -4 -6 
3.88 PESO) GED) 3.65X10 OF; OLS eS 
-4 6 
pee) bb 27X10 5% 65X10 0.029 pane | 
25 a7 
Ce Se a4 7X10 72 05X10) OFOL5 cb gee) 
5 a) 
2.00 2.06X10 im Q3X%10 0.034 halt 


Values used for thermal conductivity* 


platinum K.=,0.17447.cal/sec.cm.deg 


I 


tungsten K 0.4254 cal/sec cm deg 


*Values taken from Handbook of Chemistry and Physics, 
Boncdeedition wine: Chemical Rubber! Co. (1971) 
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expected error in) avand p, the regression) coefficients, sand stnererore 

the results of the entire study. An expected error of 8.5% was obtained. 
The average error of the Sh versus C curve was found to be 16.5%, approx- 
imately double thesexpected error. [his indicates tiatetnemtacit 1h ad 
notcontrol one of ‘the variables as well as expected, 

It is believed that the control of humidity was not as good as 
anticipated. Large errors could easily be introduced when the measure- 
ment of the humidity depended on wet and dry bulb readings. An error 
in reading “either one~or both of the temperatures couldshayve veasi ly 
created the errors shown. Another possible error source in measuring 
and controlling the humidity was the distance, between the Sourcesof 
humidity, the actual test ection, and ‘the poanuiotmneasurement..s The 
source of humidity was located approximately five feet upstream of the 
position of the supporting wire. The psychometer was located another 
foot downstreamvof this#wire. eihererore;, “ther arrimust» pass sro feet 
before its humidity can be measured. When the humidity exceeded the 
desired value and ‘the controlling circuitry tured the water atomizer 
off, an additional’ five fteet of thethumidvair passed throught themtest 
section. Additionally, the air may not have been adequately mixed and 


thus a non-homogeneous air stream may have resulted. 


Chap tere 
CONCLUSIONS AND REMARKS 

i) The experimental facility that was developed accurately and 
reproducibly measured the droplet sizes and adjusted the wind tunnel 
velocity to the corresponding terminal velocity. The heat gain from 
the supporting wire was negligible when compared to the heat transfer 
Ooutiecacropletaur hereroremathesfaci lity (simudatedstrecotalliconditions 
and represents a viable and flexible resource ,which can be applied to 
FucuTe testes Of this nature’. 

2) The procedure used to approximate the 50 micron tests was 
supported by observations of experimental results. Subsequent analysis 
of the data showed that the assumptions used in the procedure remained 
consistent throughout the study. 

3) A correlation between the Sherwood number and Renee colle 
was determined for water as: 

shietl.7551.4e80.535 Reaanse: ve 
A correlation coefficient of 0.799 was obtained with this equation and 
indicated a good fit of the data. This equation agreed well with previous 
works with the exception of the intercept term. This value was determined 
to be less than the theoretical limit and indicated a discrepancy between 
the theory and experimental -data. 

4) Percent loss of mass versus time data were obtained for each 
solution, A standard condition of 20°C, 20% relative humidity, and 
Pe Uuemicrone initia lediameter was used. Fach of the three test variables, 


initial droplet diameter, temperature and relative humidity, were changed 
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and the resulting curves plotted to show the relative effects on mass 
transfert. The sresultsmwere representat1veuor stnevwectiais proce sor ang 
an error iof—"s.5%écan be expected.” Intgenéralyetne larcesr ettrect on 
the transierrof mass sresulted irom ar changemim tnesinim 1a! dropiet 
diameter. 

S) The psychrometer measurement and the control of the air humidit; 
Wasa the “solircesot, greatest "error Min= thi s=stucyr wr more accuraceemode Of 
measurement should be employed, such as a set of wet and dry bulb 
thermistors. The source and measurement of humidity should be set closer 
to" the*position, of “thes Supporting wire with perhaps a need for oetter 
mixing. This will allow for the reduction of varying humidities along 


the wind tunnel. 


(© 
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ELIMINATION OF LOW VOLATILITY SOLUTIONS 
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Appendix 2 


COMPARISON OF SOLUTIONS 13 AND 14 TO WATER 
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Table 713 


Free Evaporation from Beakers 


Approximately 10 cm of each solution was measured 
and placed into 100 ml beakers 


Weight of Beaker and Solution, Weight Loss 


Time and Elapsed SO mee SOun LS 5017, 14 | 
Day of Time Distilled SerOZ aN ley, 1DFG2e Nels \| 
Measurement (hours) Water 100 gal. water| 100 gal. water 
4750) pm. Beaker + Beaker + Beaker + 
Sol s=59 75959 So 1 #5 7%. 95609 Sot .259.35768¢ 
1/25/32 Massmlossm= MassmlLosse= Mass Loss = | 
0.0g 0.0g 0. 0g 
8220 a.m. Beaker + Beaker + Beaker + 
15.83 Sobs=S 725522 SOli=55-59469 Sol257 751419 
82 
el Mass Loss = Mass Loss = Mass Loss = 
2.0251g 2.0071g 128627 2 


4:30 p.m. Beaker + Beaker + Beaker + 
Sol .=56.1218g Sol.=54.541g SO 1 [=5529229 


| 
1/26/82 
A Mass Loss = Mass Loss = Massmuossm= | 





3.4765¢g 3.415 3.45572 





















9:00 aim. Beaker + Beaker + Beaker + 
4 So1.=53.7291¢ 901 252755352 Sol .=53.582g 
1/277 82 
Mass Loss = Mass Loss = Mass Loss = 
S.8692g 5.602 8g S.7948g 
Beaker + Beaker + 
4 


Sol. =S1. 1825 S0l.=52 .5359¢ 










Mass Loss = Mass Loss = 
Ont oeg 6.8409g 
S=20F asm. Beaker + Beaker + Beaker + 











65585 So1.=50.2656g Sol .=48.907g 


1/28/82 





Sol .=50.2748¢ | 





MassuLoss. = Massa Lossi= Massmitossm= 
Woe S 9.049¢g 9.102g 


4:30 p.m. Beaker + Beaker + Beaker + 
0 Sol.=49.6973g So1.=48.1028g Sol.=49.4117¢ | 
1/28/82 
Mass Loss = Mass Loss = Mass Loss = 
9.901g Fess o2n 9390512 





Ulro 
4:30 p.m. Beaker + 
8.0 50.1. 2527.5154¢9 
1127782 
Ee Masse losss= 
7.0849g 
ipa 
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Table 14 


Percentage Mass Loss for 


Free Evaporation from 100 ml Beakers 


Lom@alenassm Lost 


solution IZ, = 9.9010". 
SO0LUb LONG 155 = 9853522 
SeluCi one 4) sao, 905 1. 38 


Percent Mass Loss 


Elapsed oh cee Es Sh ela ae 
Time Ssolutpon 12 oolution 13 Solution 14 
(hours) Distilled Water 5 oz NT/100 gal HO 10 oz NT/100 gal H0 
0,0 0% 0% 0% 
$5283 20,454 20737% 18.69% 
24e 352114% 34 .66% 34.68% 
40.5 59.28% 56.86% 58.15% 
48.0 71.54% 68.74% 68.65% 
O5.80 94.26% 91.84% 91.34% 


125.0 100.0% 100.0% 100.0% 
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Appendix 3 


PROCEDURE FOR APPROXIMATION OF DATA 
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I. Parameters Used in the Approximation Scheme 


D 


0 


- Original droplet diameter. Diameter of the droplet at the 


Start of the €vaporation rate history orvac times 50m into 
Value 1S) éxpressed inamicrois, 

Critical droplet diameter, Diameter at which the rate or 
evaporation of the droplet changes drastically. At this@dia- 
meter crystallizationeof the solid matter of the droplemssolurion 
begins causing the evaporation rate to decrease significantly. 
Criticalstime. “Timesiecicea tOrsthe«urup. ctulo evaporate sf com 
Do to De: Expressed in the units ef seconds. 

Final time. Time at which the evaporation can be considered 
over.” ~Expressed in seconds: 

Initial volume concentration of non-volatile material of the 


solution, a dimensionless value. 


bee 
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Assumptions Used in the Approximation Scheme 


bo 


A water-based solution evaporates at the same rate as water 


until the: solid) matierzal/) begins toricrystallaze ton:isolidify. 


tne droplet, diameter at which crystallization starts 1s 
unique for any initial droplet diameter. This value will 


not change for a variety of temperatures and relative 


humidities. 


(iesratio of the volumes of the droplet at €rystallization to 


the initial droplet is the same for any diameter droplet 


aicersthe point of “crystallizaeion, the rate of evaporation 


is not dependent upon droplet size and can be approximated by 


a’ constant) slope Jine. 
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III, Calculations for Assumptions 


Ve - Volume concentration of non-volatile material 


y. =~ Volume of non-volatile material of initial mixture 
€ Total volume of initial mixture 


Assuming each solution is completely mixed, V 
sizes, ot dropiees 


C would be constant for all 


D. - final droplet diameter, all volatile material has evaporated 


then z 
D 
4/3 ( fees) 
2 
bl AS od 
De 2 Ye 
4/3x( 57) 
or 3 
) ? Rowe 
Dy G 
SOV that 
D- 
es es = CONS ti =k A3.3 
Do G 1 


where KY is unitless. 


Using assumption 4, dD/dt - const. for t > ve for a given environmental 


condition, so that 








Dea e 
sa ieee const. = K, where dt = th chil a A3.4 
or 
Dee DAs eke ds} 
~° C 2 A3.5 
whe re K, has the units of microns per sccond. 
Combining equations A3.3 and A3.S: 
Ds 
Wer 2) ] ae A3.6 
Dy 1 
or 
D K,(dt) 
D 5 = Ky ate Sie A3.7 


13) 


Assumption 4 states that the second stage crystallization is not dependent 
UpoOmes! ce etieretorerateicanoU dependent on Size. “Examination of equation 3 


Tevedls CVeTyYChing Olathe. tightnand=sade-of-equation 5°15. constant for a 
given environmental condition. Therefore, 


D 


C A3.8 
Dp. = const. 
0 


for a given environmental condition and all sizes. This is what was stated 
in assumption 3. 


Examination of tests with different temperatures and relative humidities 


Ne reaces that the value ofredt Varics. Little betwecnatcests. Therefore, 


D 


en = COnst = A3.9 
Dy 


fOr allecmvironmentai conditions. [his 15 what was stated in assumption 2. 


Hoe 


haben. 
Summary of Errors Obtained 
with Approximation of Tests 
Statistical Formulas 
Standard deviation = o 


mean absolute error = ov2n 


probable error = 10.6745 % 


Percent 

Mean Mean Percent 

Solution Average Standard Absolute Absoiute Probable Probable 

Number D/Do Deviation EPTOY, Erros Error Error 

2 0.564 0.054 On04St 7.64% 0.0364 6.45% 

5 0.841 0.040 OF0519 3.79% 0.0269 5.21% 

6 0.400 0.018 0.0144 3.60% UFO L212 3.03% 

Z C526 O,UZ50 0.0204 3.88% O20 1S S220% 

8 OF295 OS02Z352 0.0185 6.22% OZ0156 erat: 

9 0.265 0.01488 0.01187 4.48% 0.01004 3.79% 
10 0.524 0.0519 0.0414 7.90% 0.03501 6.68% 


1M 0.465 Urvsls 0.0409 8.80% 0.0346 7.44% 


15 


Taple 16 


Parameters Required for Second Stage of Evaporation 


Average 
Solution Average second stage slope 
Number D /Dy microns/sec 
2 0.548 0.076 
: Orne 0.208 
5 0.394 0.200 
? 0.490 ee 
8 0.260 0.067 
y Doe oO 0.100 
10 0.488 OO 


unt 0.440 0.088 
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Experimental Error Analysis 


Four variables that contribute to error: 
Temperature, Humidity, Duameter, Velocity 
Errors in Measurement of the Variables 


: wi O 
1) Temperature - for a given environmental condition *.2 °C at 20eC 


sO that) TEMP = 0-01 
TEMP 





Absolute Temperature NABT = aun ees UL0GS 
ABT 293-01 








i) Humidity “HUMID = 0.05 


HUMID 


i} Diameter measurement from microscope - camera set up approximately 


= 


a : ~ 
1 micron 2 100 micron diameter ake l Se ORO 


‘ie 10 


ij Velocity - average percent error approximately 10% 


Se 0 


“SNe ikviscosityyidiffusion coefficient, and sdturation, vapor coefficient 


Vary Janerally with the air \temperature 


J 
2 s0 that “y = 0.00068 = 0.00034 


u 2 


Visc with (ABTEMP) 


= 


Del with TEMP so that V =O). 01 
V 


Satvap with TEMP so that ue = 0.01 


Nn 
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Density 
P . M 
Rete ATM a \. WwW W= 1+ 0.622 HUMID : Satva 
R + ABTEMP Parm - HUMID « Satvap 
W W W 2 a 1/2 
W  _|{ “HUMID s ik 2 2 : 
We +(e) (=) = | 9-08 + 0.01 = 0.05099 A4.1 
Ww w : W - 8 Lyd 
OL on ABT W a Z a Pr 7 
z = fe (Gj = 0.0068 + 0.05099 = 0.050995 
Error in sears 
p _ {Satvao(T 7 HUMID - Satvap(T)) 
coeff ABT, ABT, 
Ws w \2 W 2 W Zz L 
Coct ime =f CM es) Sy ff aisles ; 
Pad ece S - ABT HUMID Ad? 
ep 
i 
sweet | 2° GF 01) 7 eo (a0n UD0b8) Fe (0.05)? | : 
coeff 
WwW 
nee a 


p 02052 
coeff 





Error in Determina taoncofsSc t= me 
p 2 Ze 4 
ee (2) ea ( 2) 
eee == + cae sy 
SKe vl 0 V 
A4.3 
W 
So) (16).65197 
Se 
Error in Determination of Re on 
2 2 2 2 1. 
W Re Wy “a \ Wi ye 
——_— = — + — +{ —— +f —-— 
Re V d } fe) a 
A4.4 
Ww 
Re, - 
ee GaL2e 
Ercor in Determinations=or €. = Re ge 
W Ww 2 W r 4 
cL] VISIO Ve es 2 Yeeeah 
C {2 Re ae Na SS 
W 
aa 0.06205 
: mb rd _  -47VaM Sy - eae 
Enror ‘in Determination, of Rate, = as thi sinec ( ABT if ABT 
W : W 2 W é 4 
“RATE d v Pcoeff 
RATE etete\ tdel] | thal Vo Ay Ee 
coeff ee 


)) 
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Error in Determination of (3) = CONV. 
Cc 
2 
CONV. = onD d(D) 
2 dt 
Z 2 ue 1/2 
W W W 
2 CONSE (; ay) | a) eC 002) aeeenneon OLT 
CONV. D dD Aes 
W 
CONV. _ 
CONV oo 0.0224 
: CONV. 
= pl cathe 
EProuelieot - RATE 
W W Be Vit ate 1/2 
sh CONV. + RATE 2 2 a 
A4.8 
Shae 
arn 0.0584 


Error in 8 and a, the regression coefficients 


3 RZ 


2 
W W W 
— . (‘) f =) | = 020852 A4.9 


The error in plots of percent loss of mass versus time should be the 





Same as the error involved with the determination of the regression 
COGEii cients. 


Therefore Miler 





= 0.0852 A4.10 
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Errors Caused by the Deformation of Droplets 


b a 
Deformation of the droplets defined by the aspect ratio n=2 eh 
the volume of sphere with unit radius V = 4/3 4 Aad. 
e 2 ay Sd 
thesvolume ot spheriod V= 4/3 m7 apby = 4/3 a nye a A4.12 
equate A4.11 and A4.12: 
zs 2a 2 
Vis Vw A See a4 1 in a A4.13 
3 -2 
a = =k, weepe a= 7 oe A4.14 
replace a into the definition of the vaspect ratio 
17/3 
b=ane=n / A4.15 
find the cross sectional area 
* = 
A “= nab = oT nh ape A4.16 
equate with equivalent area of amcircle 
saa ene A4.17 
eq 
-1 
OLE Th if A4.18 
eq 
a unit radius was assumed in A4.11 
-1/6 
oie. ree =e {] A4.19 


Theretore, the ratio of the equivalent radius as determined by the 
cross sectional area measurement to the actual radius of the sphere is 
equal) to the negative one sixth power of the aspect ratio. An aspect 
ratio Up eroul.5oecan Derusedyand: still) cause only a 10% error in the 


measurement of the radius. 


Eprect Ot the Aspect Ratio on the Drag Coefficient 


Theedrag coectricicent as defined as: 


EE 2 i 
toe tee 9) y bad 
Cy (DRAG) / (1 ae” (0) V. ij A4.20 


using A4.19 Cy ae can be calculated 


aloe g o_V ) Nas 2d 


act ast 


Cy Ss = Z(DRAG)/ (7 7 
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using A4.20 and A4.21 the effect on the drag coefficient can be 


determined as: 


C 
ie WEEK Bees 44, 


Cy BIGIE 


Effece of the Aspect Ratio on the Terminal Velocity 


The terminal velocity is a function of the drag coefficient 


r eH? 
Nerd ier 


substituting A4.22 into A4.23 allows the effect on the 


terminal velocity as: 


V 
eee SEY hee 


The ‘effect ofsthe aspetteratio on the terminalavelocity, varies 
With the one sixth power of theyaspectaratios eilhis: effect. wild 


be only 10% for an aspect “ratio of'n = 0255S. 


) A4. 
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Magnetic Tape Procedures 
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The test numbers of all 496 files stored on magnetic tape can be seen 
in Table 17 along with their corresponding test variables. This table also 
lists the order’ in which these tests were stored. “The first Jineror each 
file contains six numbers that describe the test. The first of these numbers 
is the test number or name given to that file. The second number indicated 
the solution index number. The third and forth numbers give the environmental 
fyaitesesy i, in re. and relative humidity, in %, used for that specific test. 
The delay time, in seconds, for the test is displayed in the fifth number. 
Finally, the remaining number lines or records of the test is given by the 
sixth number. Each of these remaining lines contains three seconds of data. 
For each second of data, three values are given; the elapsed time in seconds, 
the droplet diameter in microns, and the actual wind tunnel velocity in centi- 
meters pen second,,, Thereforeseach line of data contains nine values, three 
times, three droplet diameters, and three velocities. 

The file numbers are determined by the testing variables of solution 
number, temperature, humidity, and initial droplet diameter. A test iteration 
number was also used to keep track of the tests when several runs were per- 
formed. The test iteration was the first digit of the test number. The next 
two digits correspond to the solution index number. The third digit represents 
the humidity used in the study. The fifth digit represents a temperature value. 
The sixth and final digit designates an initial droplet diameter. The index 
numbers and values for humidity, temperature, and initial diameter are shown 
in Table 18. This system allows a unique number to be given for any specific 
Lest. 

In order to access the files stored on the magnetic tape, a set of pro- 


cedures must be performed. These procedures along with a sample file can be 


SL 


seen in Figure 59. After logging onto the system, the main program must 
be called up. This is done by the following command; 

asg, a eVapsprogs. 
A request to mount the tape is then summitted by; 

ase. f) eVapstape., u Syv//////a;. £32032 
where £32032 is the tape number. At this point the system must be readied to 
read from tape 8 and write on tape 9. The following procedures prepare 
the system for this operation. 

uses SA evap*tape. 

asg,t evap*disk. 

WSe7o 7 evap eds K 
In order to access a given file, the system must move on the tape to the 
Start of that file. The operator may move to the desired file by using the 
command ; 

move. 8. yaN 
wheres Noisethesnumber,Or files to be moved. In the example given in Figure 59 
ben (egevOreeuhat tnesi eal) 154s stored in the twenty first position. 
Therefore, to move’ to a desired file, the file number, given in Table 17, must 
be known. Subtracting one from this file number would give the number of files 
to be moved. Once the start of a necessary file has been found, a program to 
read and write these is initiated; 

xqt evap*progs.read 
Pe istine Of this program can be séen in) Figure 60... The number of records 
contained in the file will be displayed once the program has ended. At this 
point, a temporary file has been created that will last until the operator logs 
off. In order to create a permanent file, a file name must be assigned after 


the tape has been mounted; 


is 4 


asg,up filename. 
Once the desired file has been read in, the system is located at the 
Start of the next file. If it is desired to examine an additional file, the 
operator can move from this position. In the example given, the twenty-first 
file was read in. Once this was performed, the system was located at the 
start of the twenty-second file. If, for example, the operator wished to 
examine the twenty-forth file, the move command must be used with N=2. If 
the operator wanted to examine the sixteenth file, the system may be moved 
backward with the following command; 
move,b 8.,N 
where N would be six. Moving backward tends to stretch the tape so that the 
operator may wish to rewind the tape and then move forward the desired amount 
Of filess 
The example shown in Figure S59 shows the format in which the data are 
stored in the files. The first line contains the test number, solution number, 
temperature, relative humidity, delay time, and remaining lines of data respect- 
ively. The” lines®of data consist.of three set suorselapsed Cime,=cropietraia— 


meter, and wind tunnel velocity in that order. 
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Evaporation Rate Tests Stored on Magnetic Tape 
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Pap bo eas, SOLUTION Vee HUMID INIT DIAM 
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aks Py ae UA e SOLUTION TEND HUM ED INIT OTAM 
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FILS hes Para tle SOLUTION TEMP HUMIQO INIT DIAM 
oes Tes14 1 5 90 400 
234 ay Aa 7 16 9 | 400 
255 19733% i P5 90 400 
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Te ST Mate SOLUTION TEMP HUMID INIT OIAM 
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TEST MAME SOLUTION TEMP KUMTO INIT DIAM 
aries 11 5 19 100 
Vise 1i 10 30 100 
eile ass 11 15 30 100 
US 11 20 20 100 
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SOLUTION 
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Table 18 


Index Number for Test Variables 


Index Number Relative Humidity 
1 20% 
2 60% 
5 90% 
Index Number Temperature 
1 5°C 
2 10°C 
3 15°C 
4 20°C 
Index Number Initial Droplet Diameter, Microns 
l 50 
2 100 
3 ; : 250 


4 400 
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Figure 60 


Listing of Program Used to Read and Write the Stored: Files 


EVAP *#PROGSC1).READC 7) 


CHARACTER*420 LINEC 600) 

a CHARACTERSO HDR 

o 40 CALL NTRANS (8,2, 45, HDR,14,22) 
4 Te Gee gee OO Ce Oh) 

s DECODE C€410G@@,HDRI TREC 

4 PRINT %*, TREC 

7 WRITE 097 42003 IDR 

3 AQOVQ FORMAT €50X,14,6X) 

vA naxx = RE Ge 3G 

40 CALL NTRANS €8,2,MAX,LINE,L2, 22. 
An Oe =a TREC 

42 Dein (97 44M) ENE CI) 

Sic A4M0Q FORMAT CA420) 

44 4290 FORMAT CA6Q) 

AS 20 CONTINUE 

44 PRINT naer LEE CONTAINS 21 REC, RECORDS 
47, COOSENDRALES 7 

43 St GP 

a7 END 
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